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Purpose 


The development of reliable public policies for guiding the 
growth of shopping centres and for evaluating retail proposals is 
one of most formidable and fundamental challenges confronting 
planning boards and municipalities today. Unsatisfactory municipal 
guidelines frustrate shopping centre developers and retailers alike. 


The Ministry of Treasury, Economics and Intergovernmental 
Affairs is responsible for policy aspects of The Planning Act and, 
as such, has a vested interest in the kinds of tools and techniques 
that are being used by municipalities to formulate local planning 
policy. 


This is the second publication prepared for the ministry to 
help public and private decision makers develop methods for guiding 
the planning and evaluation of major retail facilities. 


The brisk demand for the earlier report, Shopping Centre 
Dectstons: Evaluatton Guides demonstrates the clear and present 
need for such guidance. Shopping Centre Dectstons, first published 
in 1971 and since reprinted, illustrated the use of trade area 
analysis -- a step-by-step evaluation guide to decision making -- by 
detailing the results of two case studies. It pointed out 
advantages and limitations of this technique, and it urged that 
the potential of computer-based retail gravity models for use in 
Ontario be seriously explored. 


This document examines two retail gravity models which are 
capable of handling and testing information to an extent not possible 
in trade area analysis. Again, two case studies are used to describe 
the application of these measurement tools. 


Although the actual case studies were completed in 1972, 
considerable time and effort were spent on the preparation of the 
final manuscript to achieve clarity when discussing a very 
complicated subject. Meanwhile, however, there has been a free 
and useful exchange of the findings with several municipalities 
including the regional municipalities of Niagara and Waterloo, 


Oxford County and the cities of Kingston and London. 


Now that the report is available, it is hoped that any future 
applications of similar evaluations will also be published as this 
will contribute greatly to the body of knowledge of these models 
and advance their usefulness to the public and private sectors. 
Acknowledgement of the use of this publication as a source document 


would be appreciated. 


Both publications, Shopptng Centre Dectstons and Planning 
Applications of Retail Models, are the work of Robert W. McCabe, 
CD. BALM eGomeiMesc GP) a MOI Pisani McCabe is a professional 
planner with many years of experience as a consultant, university 
lecturer, municipal employee and researcher, with special 


emphasis on retail location. 


His task in preparing this document was made doubly difficult 
as he addresses two audiences -- those with experience in applying 


computer models and those without. 


Mr. McCabe was assisted by C. A. Rabenda, BA, M Sols) iemecane 
G. W. Wright, BA, M Sc(P1), who were employed by the ministry and, 
at the time, were graduate students at the University of Toronto's 


Department of Urban and Regional Pilani un. 


Acknowledgement is gratefully made of the cooperation of the 
many municipal people who helped in the preparation of material. 


Mort wesainclaiy, eNerr 


Manager 
; Special Studies Section 
Toronto, 1974 Local Planning Policy Branch 
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1 Introduction: Bases, Approaches and Conclusion 


The primary purpose of this study has been to explore the 
potential of gravity type models to assess the impact of new or 
expanded suburban shopping centres on the existing retail commercial 
facilities of a municipality. The study is a sequel to an earlier 
Study published in 1971 entitled Shopping Centre Decisions: 
Evaluation Gutdes (McCabe 1971). 


The previous study dwelt at length on the elements of risk 
inherent in shopping centre decisions and the adequacy of trade 
area analysis to answer the questions that decision makers must 
resolve to their satisfaction. 


This study has been aimed at amplifying what is known about 
retail models and their ability to answer these and other related 
questions with equal or even greater reliability than current 
methods allow. 


Retail Models and Shopping Centre Decisions 


The potential usefulness of retail models for urban and 
regional planners and municipal decision makers rests on the claim 
that models can be helpful in providing answers to such basic 
planning questions as: How much retail space is required? Where 
Should it be located? When will it be needed? 


In more detail, retail models attempt to answer the following 
questions about shopping centre decisions: 


1 What potential is there for existing or new shopping centres 


in the area? How many centres and what size of centre can be 
supported? Which locations are best? When will new centres 
be required or when should existing centres be expanded? 


How much in annual sales will be gained or lost by each centre 
when a new centre is introduced or an existing centre expanded? 


Which pattern of centre sizes and locations offers the greatest 
potential or has the least adverse effects? Which pattern is 
most likely to maximize overall customer shopping satisfaction? 
Which best provides comparable levels of service to all 


residents in the area? 


What are the consequences of alterations in the assumptions 
made about future travel times and distances, levels of per 
capita income, population growth rates, distributions and 
densities, and how might each affect retail location and 


shopping demand? 


How similiar is the set of future centres to the set that might 
result from normal market processes? How Meafficient is sche 
planned system of centres likely to be in producing reasonable 
productivity for the retailer measured in terms of centre size 
and sales per square foot, and how effective is the system 
likely to be in minimizing average trip lengths for the 
shopping public? 


What levels of consumer flow are most likely from each 
residential area or neighbourhood to downtown and suburban 
shopping centres? What levels of sales are generated at each 
centre? How closely do these approximate existing conditions? 


Like any method which attempts to evaluate what is likely to 


happen jn the future under varying circumstances, retail models 


have their limitations. In the case of trade area analysis, which 


appears to be the only currently available alternative to retail 


models, the principal limitations are: 


ie 


the high cost in time and money needed to test even a limited 
number of alternative arrangements under a restricted variety 
of assumptions about the future, 


the lack of precision concerning the uncertainties involved in 
the use of estimates, and 


the large number of factors that have to be kept straight in an 


analyst's mind while he is doing his step-by-step calculations. 


Retail models appear to overcome these problems by 1) allowing 
for the evaluation of an almost unlimited number of possibilities 
at low cost, 2) by delimiting precisely the effects of weak 
estimates or poor judgement, and 3) by making use of the ability of 
the model to keep large numbers of factors in focus while the 
necessary calculations are being made. 


Retail Models and Public Policy 


Retail modelling also appears to be the best means currently 
available by which municipal, regional and provincial governments 
can arrive at answers to policy questions crucial to their ability 
to plan for orderly growth and development. 


Official plans in Ontario must by their nature incorporate 
provisions to ensure that residents will continue to enjoy the 
benefits of reasonably located shopping facilities. 


Retailers too must be assured of a consistent set of policies 
within which they can plan the best means of serving their customers. 


Failure to face these important policy issues invites the 
possibility of less than satisfactory conditions for shopping 
centres' consumers, retailers and property owners, and creates 
serious problems for elected representatives and their officials. 


Among the major policy questions that retail gravity models 
can help answer are: 


ds What are the probable consequences of a lack of policies or of a 
continuation of current policies and practices on the location 
Oia nretailefaciiicies? 


ae Are policies on retailing consistent with other policies 
concerned with population distribution, densities, and the 
provision of transportation and piped services? 


3% To what degree do existing retail facilities adequately serve 
existing populations and what will be the effect of trends to 
new forms of retailing? 


4, At what scale do retail facilities cease to be of purely 
municipal concern and become a matter of regional interest? 


a 


5 To what degree does a proposed shopping centre, or several 
proposed centres, match the future desirable state of affairs 


envisaged by the official plan? 


6. How much flexibility in defining exact retail locations can 
safely be built into a plan for the future? 


Pressures for new and better located retail facilities will no 
doubt continue to force elected bodies at municipal, regional and 
provincial levels to seek ways of answering these and other policy 


questions. 
APPROACH AND OUTLINE 


The approach adopted for the evaluation of retail models is 
based on the same assumption that was made in Shopping Centre 


Dectstions: Evaluation Guides. 


Decision making is a matter of selecting a solution to a 
problem from among a number of perceived alternatives. The adopted 
solution is hopefully one that involves the minimum risk for the 
decision taker. 


The planning process is designed to help the decision taker 
-- entrepreneur or politician -- make sure that risks are minimized 
by removing as far as possible uncertainties about the future. 


This study has revolved around the ability of retail models to 
reduce uncertainty and to evaluate a wide range of possible 
alternatives under varying assumptions about the future in a 
realistic and pragmatic manner. 


Chapter 2 reviews the development of retail models of the 
gravity type and describes in some detail the two models selected 
for intensive study. Understanding these models and how they 
operate should help the reader comprehend the results of the 
analyses contained in later chapters. 


Some time is spent considering the conceptual, operational and 
informational bases of the two models so that a detailed comparison 
may be made between the legitimacy of retail models and of the only 
other means of evaluation, trade area analysis. 


It has been assumed that readers are familiar with the way in 
which provision for orderly growth and development is achieved 


aeha 


under The Planning Act in Ontario and the role that has been played 
by market studies in the past. 


Chapter 3 presents two case-study applications in Ontario 
employing the two models selected for examination. 


Concern over the growing number of applications being received 
by local municipalities for permission to construct suburban 
Shopping facilities appears to be most critical in urban centres 
with populations from 25,000 to 150,000. 


Studies of the changing structure of retailing in North America 
indicate that the prime tenants in the larger shopping centres are 
currently moving into smaller communities in a major building 
expansion. Both case-study municipalities typify this situation. 


Large centres on the other hand generally have access to the 
research and informational resources needed to carry out major retail 
studies. At the same time the addition of another centre in a large 
metropolitan region is less likely to have the critical impact of a 
large suburban shopping facility in a smaller community. 


Finally, the complexities of retailing in very large centres 
such as Metropolitan Toronto might obscure the fundamental lessons 
to be learned about retail models. 


A testing procedure is developed here for answering the 
questions posed in Chapter 1. Conclusions are reached by comparing 
the two versions of the retail model tested and the results achieved 
in each of the two case-study municipalities. 


Chapter 4 contains a general review and evaluation of the 
results of the study. Special attention is given to aspects of 
retail modelling that require particular care by those who provide 
the information on which the models operate. The basic information 
used in retail models is shown to be largely identical with that 
required for trade area analysis. 


A Technical Appendix has been prepared as a companion piece 
to the report. It sets out the operating instructions for the two 
computer models employed in the study in enough detail to permit 
their subsequent application at a municipal or regional government 
level by planning staffs or consultants. 


Since computer language and technology mean little to the layman, 


the appendix attempts to explain in non-technical language the 
procedures used in the two programs -- particularly for those who 
must place high reliance on the results of computer models. 


Fach of the two programs is presented separately. The 
information that is entered into the models (the input variables) 
is described and the information in the printed results is listed 


in detail (the output variables). 


Program listings, explanations of the exact nature of each 
step and a flow diagram are provided for both programs. Additional 
options that the user may wish to use under certain circumstances 
are also explained and the necessary modifications in control cards 


or program decks specified. 


The Annotated Selected Bibliography contains those publications 
on retail models -- many of them technical -- that users would 
likely want to become familiar with before attempting to apply 
either of the retail models or to interpret the results. 


The General Bibliogzapny lists alll books, journal articles and 
government documents that were consulted or referenced during the 


study. 
SOME GENERAL CONCLUSIONS 


The application of retail modelling in the two case-study 
municipalities appears to have substantiated the powerful analytical 
advantages claimed for the models tested. 


This was particularly true when they were employed in 
conjunction with conventional step-by-step trade area analysis. In 
both instances and for both models a close replication of existing 
retail shopping patterns was achieved. The dependability of short- 
term forecasts made by the models was amply demonstrated. 


The models provided reliable and usable answers to the planning 
questions that were posed at the outset of the study and to the 
policy questions. 


Differences were found, however, in the degree of reliability 
and the ease of operation of the models employed. 


The Huff model proved to be easier to put into operation but 
less satisfactory in achieving an immediate close replication of 


actual conditions. 


The Cullen model proved to be superior in the manner in which 
the results are presented, in the reliability of its forecasts, and 
in the ease with which a close approximation of observed results is 
obtained. 


Extensive applications of the models in as many operational 
situations as possible would now seem appropriate. It is important, 
however, that the results of all applications be published as a 
matter of course so that all those interested in the development of 
public policy may benefit. This is the only way expertise can be 
developed and the benefits of more widespread use realized. A 
number of these studies have already been made and more are underway 
at the time of writing. 


More detailed conclusions will be found at the end of each 
chapter. 


EMERGING ISSUES 


Discussions with planning agencies and others about the nature 
of the studies undertaken brought out some further problems 
requiring attention: 


We Where regional governments exist, studies need to be undertaken 
at a regional scale even though a proposed centre may be 
located in a single constituent municipality. 


2. Where regional governments do not exist, some co-operative 
means must be found to allow broader analysis of the 
implications of establishing suburban shopping centres in 
those municipalities that lie immediately outside the large 
urban centres. 


B. An assessment of the implications of new forms of large-scale 
retailing that are just now emerging may be difficult to 
determine using only the currently available methods. Further 
refinements of trade area and retail modelling techniques may 
well be needed. 


It would seem logical that the lead in developing new approaches 
to evaluating retail locations should come from the Ministry 
of Treasury, Economics and Intergovernmental Affairs, either 
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through its own staff or through arrangements with local or 


regional planning agencies. 


2 Retail Modelling 


Two formulations of the retail model have been investigated in 
the course of the research undertaken in this study. 


The first is the Huff and Blue version of the Huff model which 
has been extensively revised in the Department of Urban and Regional 
Planning at the University of Toronto and, later, by the study 
researchers. 


The original program is fully described in A Programmed 
Solution to Estimating Retail Sales Potentials (Huff and Blue 1964). 
The revisions were introduced to improve the tables produced by the 
model and make them more useful for planners and comparable to 
those produced by the second model. 


The second formulation is the Cullen version of the Wilsonian 
entropic model (Cullen 1969) also extensively revised at the 
University of Toronto by the study researchers. A comprehensive 
understanding of this model requires a reading of the papers by 
Hayes (1968), Wilson (1970b), Batty (1972a) and Cullen (1969). 


This chapter is written for those who want to gain a general 
understanding of the development and use of the two retail models 
used here without having to read all the relevant technical journals. 


Previous approaches to an understanding of the distribution of 
shopping facilities and consumer shopping habits are reviewed, and 
the various models that have been developed are described. A 
comparison is made between the principal characteristics of the two 
versions tested. 


This rather detailed description of the two models is 


considered a necessary preliminary to appreciating the analysis 
reported in the case study chapter. 


DESCRIPTIONS OF RETAIL SHOPPING ACTIVITIES 


Early Gravity Models 


Early investigations of the market area concept, relating the 
purpose, frequency and length of shopping trips to a hierarchy of 
centres, closely paralleled a set of empirical studies on spatial 
interaction. The underlying theme was that a definable relationship 
exists between the size of a city and the distance from which 
shoppers are attracted. 


The gravity concept of human interaction postulates that there 
is an attracting force between two areas of human activity. This 
force is created by the mass of the two populations and by a friction 
against interaction that is exerted by the intervening space over 
which the interaction must take place. Interaction between two 
centres is seen as varying as some function of population 
concentration and the distance between them. 


It was not until W. J. Reilly postulated his "Law of Retail 
Gravitation" in 1929 that the analogy with Newton's Law of 
Gravitational Force was specifically applied in marketing (Reilly 
1929). Reilly first formulated his gravity concept to describe the 
attraction that two adjacent cities exerted on intervening small 
communities for the purchase of fashion goods. Reilly's Law states 
that: 


Two centres attract trade from intermediate places 
approximately in direct proportion to the size of the 
centres and in inverse proportion to the square of the 
eta from these two centres to the intermediate 
place. 


In notational form the relationship is: 


where Us and Ui, represent the fashion goods trade drawn to two 

cities A and B from the intermediate place, Bes and Dy, the distances 
from the intermediate place to A and B, and Ee and Pb the populations 
of the two cities. 


Over the next decade the model was revised and broadened in 
its applicability. A break point was postulated between the trade 
areas of the two cities. This was the point at which customers 
were indifferent as to the city they shopped in. The formula was 
consequently modified to read 


eee ab 
Pa 

Pee 
Pb 


where dD. denotes the distance from city B to the line of equal 
attraction, Dab the distance between the two cities, and He and Ph 
the populations of the two cities (Fig. 2-1). 


This formulation, like most of the early quantitative approaches 
to marketing problems, was only one of several attempts at deriving 
exponents or indices to explain a variety of distance-related 
phenomena. 


The power function of 2, a squaring of the inverse distance, 
was arrived at by Reilly from field observations. The observed 
values, which fell between 1.5 and 2.75, were simple averages. 
Reilly made 255 field measurements in all and found that some 
exponents ranged from 0 to as high as 12.5. It has been claimed 
that for this reason alone Reilly's Law had a weak empirical basis 
-- an unvalidated hypothesis arrived at intuitively through an 
analogy with pure physics. 


Marketing Analysis 


Although Christaller is credited with hypothesizing Central 
Place Theory in 1933, formalizing the concept of a hierarchy of 
centres, the discipline of marketing had developed techniques for 
describing market activities much earlier. 


Theories of retail location evolved based on estimates of trade 
areas for goods of different categories. The concept of convenience 


= hb Re 


REILLY'S BREAK POINT ANALYSIS 


Determining the Breaking Point 


The following formula is derived from Reilly's Law of Retail 
Gravitation. It has been used for determining the boundaries of 
a trade area, the 50 percent breaking point -- the point at which 
a town secures one-half of its trade. 


At the breaking point two cities attract retail trade 
“approximately in direct proportion to the population of the two 
cities and in inverse proportion to the square of the distance from 


these two cities to the intermediate town". 


A - Larger town 
B - Smaller town 


Breaking Point, Miles from B = Miles between A and B 


hoe Population of A 
Population of B 


Trade Area of A Trade Area of B 


IZ miles 6 miles Town B 


Population 
120,000 


18 miles Population 


30,000 


Fig. 2-1: Reilly's Break Point Analysis 


Converse, Paul D. and Heugy, Harvey W., 
Elements of Markettin omcnd Edo. New York. 
Prentice-Hall Inc., T9a0, pp 641-2, 


STEP-BY-STEP APPROACH TO TRADE AREA ANALYSIS 


STEP. I Establish a tentative trade area for the 
centre 
STEP II 


SAS 
SEP lV 


STEP V 


SEP VI 
saeP VII 


Fig. 2-2: Step-By-Step Approach to Trade 
Area Analysis 


Determine the initial desirable size and 
nature of the contemplated centre to 
assess its possible power of attraction 


Estimate the total potential purchasing 
power of the residents of the delineated 
trade area for the goods and services 

represented in the centre 


Analyse the competing centres in the trade 
area and their relative attractiveness in 
relation to the proposed centre 


Survey consumer shopping habits in the 
trade area and determine the characteristics 
and attitudes of the shoppers, centre 

Datronage, cc. 


Estimate the share of the potential consumer 
expenditure in retail goods and services 
which will likely be attracted to the centre 
and translate this into retail) floor space 

requirements by type of goods 


Evaluate the findings of the trade area and 
market analyses in relation to the initial 
size and floor space allocation in the 

contemplated centre and adjust as necessary 


McCabe, Robert W. Shopping Centre Decisions: Evaluatton 
Guides, Toronto, Province of Ontario, Ministry of Treasury, 
Economics and Intergovernmental Affairs, 1971. 


buying and comparison shopping as delineators of market areas and 
centre sizes existed in marketing literature prior to 1920 (Copeland 
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Marketing studies first concentrated on problems dealing with 
individual shops or classes of goods and their market areas. 
Advances in marketing techniques continued to be made and soon whole 
centres were being analysed instead of just individual outlets. 


Later marketing studies of this nature came to be known as 
Trade Area Analysis or the Step-by-Step approach. The methods 
currently employed have remained substantially unchanged from the 
form it had developed by the 1930s. 


Trade area analysis is discussed at length in Shopping Centre 
Deetstons: Evaluation Guides (Fig. 2-2). Until very recently this 
method of analysis dominated studies of urban and regional markets. 


The assumptions underlying trade area analysis are not 
dissimilar to those of Central Place Theory and early gravity 
formulations: 


iz Primary and secondary trade areas are delineated and 
Proportions of the consumer expenditures in each allocated to 
individual shopping centres. Relatively exclusive trade areas 
are assumed approximating a closed System. 


2. The analysis is centre oriented. 


3. A hierarchy of centres is postulated based on the types of 
goods carried. 


4, The closest centre appears to be the most attractive. 
oh The size of the trade area varies with the size of the centre. 


6. Centres are evaluated one at a time. The evaluation of 
alternatives is time consuming and involves many value 
judgements on the part of the analyst based on his experience 
and. the proprietary information held. 


Trade area analysis has focused almost exclusively on the 
concept of retail trading areas like other early retail shopping 
Studies. It appears to have been Satisfactory in fulfilling the 
intended purpose of evaluating the retail potential of a single 
site. Theoretical and practical limitations appear when it is used 
for other purposes. 


ee ore 


a 


A Probabilistic Approach to Retail Modelling 


It is only recently that tentative attempts have been made to 
develop a body of theory that links retail location studies to the 
major advances that have been made in marketing in the study of 
consumer behaviour. Early methods did not reveal why observed 
regularities in shopping patterns occurred nor were the hypothesized 
levels of activity based upon an accepted theory of how consumer 
shopping choices are actually made (Luce 1959). 


Empirical studies in urbanized areas using customer shopping 
surveys indicate that consumers do not appear to discriminate among 
choices perfectly. For example, consumers do not appear to select 
one alternative consistently. The lack of consistency in choice 
becomes especially apparent where there are a large number of 
shopping centres of various sizes within a reasonable distance. AS 
a result, distance in urban areas appears to provide significantly 
less protection from competition. A shopping centre is seldom able 
to dominate an exclusive trade area. 


Trade areas also overlap because consumers find it pointless 
to continually discriminate among centres when the perceived 
differences between them are small. 


A consumer may be uncertain as to the ability of each shopping 
centre to fulfill her various shopping trip expectations. She may 
visit one or all of them at some time or another in an information- 
seeking process. 


But the next time she is confronted with the same choices in 
selecting a shopping centre, she will probably tend to choose among 
them in some more consistent manner. 


This observation of shopping behaviour, especially in a highly 
urbanized area, differs in its point of emphasis from the traditional 
trade area analysis approach which begins by examining the site of 
the shopping centre. One of the strongest features of retail model 
formulations is that they do take this characteristic of consumer 
behaviour into account explicitly. 


This new approach to the problem of dealing with trade areas 
was outlined in the early 1960s by D. L. Huff (Huff 1963). The 
apparent strength of Huff's probabilistic view of human choice 
behaviour over previous methods and the power of this version of 
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the retail model encouraged widespread research into retail models 
in the U.K., Scandinavia and elsewhere (Rhodes and Whittaker 1967; 


Agergard 1969). 


Huff suggested that the areas of competition between shopping 
centres continually overlap. Consequently, shopping activities 
appear to be apportioned among the accessible centres in a 


probabilistic manner. 


If the apportionment is probabilistic, the total of al] 
shopping probabilities must be unity. If there are four centres, 
for example, and the probability of consumers in a given location 
of shopping at centre A is 0.5, at centre B 0.3, at centre C 0.05 
and at centre D 0.15, the sum of the probabilities will equal 1.00 
as required by the theory of probability. 


DESCRIPTIVE AND NOTATIONAL FORMS OF RETAIL MODELS 


Description of the retail models is split at this point into 
two sections. Part A is purely descriptive and is provided for 
those who are not fully familiar with mathematical notation. Part 
B is in notational form for those who find that means of expression 
simpler to follow. 


A. The Descriptive Model 


(1) The probability that a customer who lives in a certain area 
of a community, defined as a population zone, will be 
attracted to any one retail centre is the ratio of the 
attractiveness of the centre to the distance from the zone 
centroid to the centre raised to a power, defined as lambda 
(A), divided by the sum of all the similar ratios for all the 
zones and all the centres. 


(2) The sum of all these probabilities must sum to unity, that is, 
WO" 


Sh 2 


B. The Model in Notational Form 


(1) The probability that a consumer from zone i will be attracted 
to retail centre j can be expressed as: 


A. 
dank 
44 
Pig z= SLE Ba (for any zone) (1) 
Seed ny 
a lee) 
where 
ij = the probability that consumers in zone i will be 
attracted to centre j; 
Aj = the attractiveness of centre j; 
ij = the distance from zone i to centre j; 
r = a parameter value which is to be estimated 


empirically to reflect the effect of travel 
distance, time or cost on shopping trips. 


fey. subject to 


n 
ap Ue er kd (for all zones) (2) 
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A. The Descriptive Model 


(3) The amount of retail sales that are available in any zone is 
the product of the number of people who reside in the zone and 
the average per capita expenditure on retail goods and services. 


(4) The portion of the available retail expenditure from the zone 
that will be attracted to any given centre is the product of 
the total available expenditures and the probability that 
customers from that zone will be attracted to the given 


centre. 


(5) The sum of the sales attracted by each centre from a zone will 
equal the total available retail expenditure in that zone since 
the total of all the probabilities must equal unity. 


If the sales at each centre from each of the zones are summed, 
i.e... the sales at Centre A from Zone |. pluses the sales 10m 
Zone 2, plus the sales from Zone 3, etc., an estimate is 
obtained of the probable sales at that centre. 


In the four-centre case mentioned earlier, if the total 
available retail expenditure is $100,000 -- the product of a 
per capita retail expenditure of $500 from 200 people -- sales 
at Centre A from that zone will be $50,000 ($100,000) x 0.5), 
at Centre B $30.,000,,at\ Centre € $5,000 and ateGentresD 
$15,000. 


(3) 


(4) 


B. The Model in Notational Form 


The total available retail expenditure from zone i can be 
expressed as: 


Eee aca (3) 
where 
E. = retail expenditure in zone i 
B. = average annual retail expenditure per capita 
C. = the number of consumers in zone i 


The amount of business attracted to a centre j from zone i 
can be expressed as: 


ar a) ere eee (for any centre) (4) 


344 = the retail sales attracted to centre j from zone i 


Substituting, it can be seen that 


A. The Descriptive Model 


If the sales of each centre are summed, that is, the sales at 
Centre A plus the sales at Centre B plus the sales at Centre C 
and so forth, the total of all the retail sales by all the 
centres is produced. 


If the total available retail expenditures for each zone are 


summed -- the expenditures in Zone 1 plus the expenditures in 
Zone 2, etc. -- the total of all retail expenditures is 
obtained. 


By definition the total of all retail expenditures by customers 
must equal the total of all retail sales made by stores. 


It is obvious that some people may spend some of their money 
for retail goods outside the community altogether. An 
adjustment must be made for this outflow. It is equally 
obvious that some people making purchases in the stores in the 
community may not be residents. An adjustment must be made for 
this inflow as well. 
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B. The Model in Notational Form 


(6) Subject to the constraint that 


n n 
yee ae (for all centres) 
Si Rs a 
where 
Sij = the sales attracted to centre j from zone i 
E. = the expenditures made in zone i 
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A. The Descriptive Model 


Having calculated the estimated share of sales coming from 
each zone and arriving at an estimate of the sales at each 
centre, it is possible to compare this estimate with the 
actual sales at each centre. Actual sales are not normally 
available from centre operators, but methods have been devised 
for estimating these independently of the steps followed 
previously. If the calculated sales and the estimated sales 
do not match within reasonable limits it is assumed that the 
difference is accounted for by the friction of distance and 
that adjustment for this factor can now enter the calculations 
by varying the value of lambda (\).9 “By thease ofe@search 
techniques it is possible to arrive quickly and simply at a 
value for lambda (\) such that the calculated and actual sales 
are in close approximation. The model is then said to be 
calibrated. 
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B. The Model 


in Notational Form 


In a hypothetical case, a summary of the site information 


produced by the model might show the following results of 
calibration at various values of lambda (i) in table form. 


Lambda (A) = 0.8 


Calculated 


Centre Sales 
($ million) 
] Z0R2 
2 7.8 
3 54.9 
4 M2 
5 8 
Totals 97.9 
Lambda (A) = 0.85 
Calculated 
Centre Sales 
($ million) 
] 20) 
2 five) 
3) 55.0 
4 9.3 
5 5.6 
Totals 97.9 


Year = 1966 
Actual Difference 
Sales Squared 
($ million) (ex al Oise) 
20.6 126 
Th ks) - 
54.8 0.1 
8 gab 
) Oz 
97.9 34 
Year = 1966 
Actual Difference 
Sales Squared 
($ million) ( x 10!1) 
2026 205 
Thats 0.1 
54.8 nA 
8.8 225 
9 0.9 
97.9 6.4 


In this case the first table appears to achieve a closer fit 


between calculated and actual sales and the calculated value 


of A = 0.8 is more probable. 
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A. The Descriptive Model 


To test the validity of the calibration, the practiGe Sto 
make a forecast or a backcast using the model as it calibrated 
to see how closely the forecast comes to conditions that are 
known to exist at the forecast or backcast date. For example, 
if the model is calibrated on 1966 information, an attempt can 
be made to predict what the model says the 1971 situation 
would be, given more population, higher incomes, new or 
expanded centres, etc. The forecast for 1971 is then compared 
with known 1971 information and the predictive accuracy of the 
model judged. If the predictions are faulty, a re-examination 
can be made of all the factors that have entered into the 
model for each year to see if there are any errors or bad 
assumptions. Generally, the predictive power of the model is 
found to be very good once inadvertent aberrations are removed. 
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B. The Model in Notational Form 


In a hypothetical case a comparison of the estimates for 1971, 
with lambda (A) set at the 1966 best value of lambda (A) = 0.8, 
might show something like: 


Lambda (A) = 0.8 Year = 1971 
Estimated Actual Difference 
Centre Sales Sales Squared 
(x 10!!) 
] 22.4 22.9 2.5 
2 9.5 9.7 0.4 
3 We ale 56.4 32.6 
4 Loe 9.8 0.9 
5 9 6.2 0.9 
6 2 Sel Om4 
Totals Teel ibe eat ae | 


OAs = 


(9) 


AveinesDescriptive Model 


A refinement to the Huff formulation, known as the Lakshmanan 
and Hansen model, was devised by Lakshmanan (Lakshmanan and 
Hansen 1965). Further extensions have been made by a variety 
of researchers. Lakshmanan made the adjustment to raise both 
the attractiveness factor and the distance factors to a power. 
In this case these are known as the alpha (a) and beta (8) 
values. The model is calibrated as before but both a and B 
values may be adjusted to bring the calculated and actual sales 
into line. (B is identical to A’ inthe earlier model. ) 


SP 


B. The Model in Notational Form 


(9) Where both A; and dij are given parameter values that enter 
into the calibration, the Lakshmanan and Hansen formulation 
reads: 


7 B 
Coen ee (7) 
1c) te eee Aad 
A re 
aS VJ 
subject to the constraint that 
n n 
eS eit See Ee (forall ecentres ) (8) 


where 


a and B are two parameters, one of which (a) reflects the 
variations in attractiveness caused by assortments, 
differences in retail shops, etc., and the other (8) 

is identical with (A) in earlier versions of the model. 
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A. The Descriptive Model 


(10) Wilson has suggested that the retail version of the general 
gravity model may be more rigorously formulated through the 
processes of statistical mechanics. In this form of the model, 
an exponential function is used in place of the power function 
of a and g employed by Huff and Lakshmanan. The version of 
the gravity model used during this current study was developed 
by Cullen at University College, London in 1969 and employs 
Wilsons! arguments. Newer developments in the retail modelling 
field typically rely on similar, more rigorous formulations. 

It was thought best, therefore, to adopt the Cullen version 
from the start. In effect, the Wilsonian approach overcomes 
any of the criticisms of the retail model as having a weak 
empirical basis. 


Earlier versions were founded on a tenuous analogy between 
retail shopping patterns and Newtonian physics. The Wilsonian 
version is derived on a strict empirical basis in a rigorous 
mathematical manner. 
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B. The Model in Notational Form 


(10) The Cullen version of the Wilsonian entropic model takes the 


form: 
SE pees (3) 
ij n : 
5 ow Adis 
Je 
where 
ij = sales attractiveness of zone i to centre j 
E. = the household or per capita retail expenditure in 
zone i 
Wy = an index measure of attractiveness of centre j 
dij = some measure of cost of travelling from zone i to 
centre j 
O. = a shopping centre attractiveness factor 
r = a travel cost parameter 
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THE SELECTION OF VARIABLES 


The ability of a shopping centre to fulfill the shopping 
demands of a consumer is obviously determined by a host of factors. 
Among those substantiated are newness, the presence of an enclosed 
mall, air-conditioning, the reputation of the retailers represented, 
the size, choice and price of goods, ease of parking, distance from 
home and advertising effectiveness. Attempts to weigh the individual 
importance of such factors have met with limited success, 


Recent research has confirmed that only a few of the variables 
are important enough to be included in workable models that ensure 
valid and reliable predictions. In food Supermarkets, for example, 
price, location, quality and service appear to dominate all other 
marketing factors (Tigert 1971). 


In the retail models developed to date, two principal variables 
are assumed to exert enough influence that collectively they can be 
safely used as proxies for all the other variables in making 
predictions. These two variables are: 


(a) the number of items that are carried by a shopping centre, and 


(b) the travel time involved in getting there. 
Items Carried 


A consumer cannot always be positive in advance that a 
particular shopping centre will fulfill the purpose of the trip 
being planned. The shopping centre may or may not have the 
particular item or brand of item that a consumer wishes to purchase. 
However, a customer usually has some prior knowledge that a specific 
shopping centre will satisfy her needs and presumably this is based 
on the number of items she feels are carried by the centre. So the 
greater the number of items carried the greater will be the 
expectation that a shopping trip will be successful. 


Consumers, therefore, seem to prefer larger centres and to 
Show a willingness to travel further to reach them. The degree of 
uncertainty in consumer shopping decisions appears to diminish as 
the size of the shopping centre increases, reflecting the assumed 
increase in the number of items available. 
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Distance 


The desirability of a shopping centre is also assumed to be 
influenced significantly by the effort and expense involved in 
getting there. The consumer bears the burden of movement to the 
centre in terms of the required preparation time, monetary costs 
and travel time invested. Consequently, there would appear to be 
good reason to suggest that a consumer attempts to minimize these 


costs in some way. 


However, there is strong evidence that the minimum cost choice 
is not consistently selected to the point where shopping is done 
only at the nearest centre. The desire to minimize the costs of 
movement reflects the notion of opportunity costs. The consumer is 
assumed to have only so much time. If the time is devoted to more 
than one activity at the expense of others, the costs of foregone 
opportunities rise to act as a check against further time losses. 
Therefore, the value of a shopping centre to a consumer seems to be 
jn some way inversely related to the effort and expense involved in 
getting to it. 


In many instances, especially for jntra-urban shopping trips, 
there are only minor differences in distance to each of the centres 
that appear to the consumer to be accessible. In these circumstances, 
the effect of distance on the selection of any one centre over the 
others would be minor. Current retail model formulations 
realistically reflect this thinking. 


Moreover, willingness to travel is not constant for all items. 
The more frequently purchased convenience goods do not command as 
great a willingness to travel on the part of the consumer as do the 
less frequently purchased, higher cost shopping goods. 


Different types of merchandise appear to require different 
amounts of consumer searching. Variations in the effect of travel 
time for different kinds of shopping trips are reflected in very 
recent retail models that deal with each level of centre 
hierarchically. Models of this type appeared only after the 
current research was completed and have yet to be released for 
evaluation and testing. 


Parameters 


The majority of early efforts by Reilly and his followers to 
delineate trade areas had assumed that a power function of 2 -- the 
Squaring of the distance -- was inherently correct and that this 
value would apply equally to intra-urban and inter-urban shopping 
trips. 


This assumption became the target for much criticism when 
studies revealed that parameter values. appeared to fluctuate 
depending on the type of trip being analyzed as well as the 
geographical area in which the study was made. 


Parameters are used in most models to equate predicted to 
observed results. A parameter is a weight or effect not directly 
measurable which may or may not be estimated independently of the 
model. Distance and attractiveness parameters are examples in the 
two models used in this study. 


Where parameter values cannot be estimated independently, they 
must be determined internally in the model by a calibration process. 
For example, where attractiveness is being measured according to 
some direct scale of shopping provision such as items carried, 
floorspace or sales, there is a case for introducing an exponent 
to allow for the extra attraction of choice. 


In some instances it has been observed that two shops of a 
given type in one centre may be more than twice as attractive as a 
single shop of the same total size, presumably because of the 
cumulative attraction offered by the added choice in shopping. The 
exponent would take a value other than unity in this case. 


Variations in parameter values may be assumed to allow for 
factors such as the added attractiveness of a large centre over and 
above its size to a smaller centre because of its wider choice of 
goods, ease of comparison, and other intangible factors not 
incorporated in the model. 


Examples of factors which are usually present in varying 
degrees for different size centres are: 


ie an apparent increase in amenity as centre size increases, 
Zt the superior attractiveness of major "name" food and 


department stores usually found in a larger centre, 
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on greater ranges; assortments , depth of stock, etc. 


Where attraction measured in terms of size alone is being used 
as a proxy variable for a number of other attributes, the use of the 
attraction parameter reflects the accumulation of these unmeasured 
factors during the calibration process. But not all models 


incorporate this particular parameter. 


A distance exponent is introduced in all model formulations to 
reflect the conclusions of earlier studies that different types of 
merchandise support differing amounts of consumer effort. The 
problems created when two-parameter versions of the model are 


employed are discussed later. 
THE CHOICE BETWEEN MARKET STUDIES AND RETAIL MODELS 


There are two approaches to the study of intra-urban or intra- 
regional retail systems of centres. One entails subjective 
judgements based on observations of similar conditions from which 
generalizations have been drawn, for example in the step-by-step 
method. The other makes use of empirically derived mathematical 


formulations, as do retail models. 


Mathematical models which incorporate probability concepts 
have a number of important advantages over subjective methods of 
evaluation for the urban and regional planner. 


lp The real world is represented in abstract and simplified terms. 


2 The extensive analytical capabilities of the computer can be 
effectively used. Large scale quantification is essential 
when trying to estimate the effect of one or more changes on 
an area containing many centres. The problem that an analyst 
faces in trying to keep a multiplicity of interacting 
variables straight in his mind at one time is overcome. 
Flexibility is thereby enhanced and a large number of options 
can be tested. 


3" The speed of computer applications produces instant feedback 
and provides a methodology that is relatively inexpensive and 
easy to use. 


4, Acceptable solutions can be reached more quickly by allowing 
the model to eliminate those proposals that do not meet 
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predetermined criteria. 


5s The large amounts of statistical data available for urban 
centres can be readily analyzed and speedily retrieved in 
readable form. 


Like the step-by-step approach, the probability models 
described in this report are static in nature. They examine only 
successive states of temporary equilibrium over time. The 
intervening processes of change are not taken wholly into account. 


NECESSARY ASSUMPTIONS IN OPERATING RETAIL MODELS 


A simplified block diagram of the retail modelling process is 
reproduced from Shopping Centre Decisions: Evaluation Guides. 
(Fig. 2-3). Five general areas of interest underly the operation 
of retail models. These are the boundaries of the study area, the 
Origins of shopper trips, measures of attractiveness and of 
distance, and locations of new centres. 


Study Area 


Unlike previous methods, probability models do not assume that 
a clearly delineated closed market area exists for any one centre 
or system of centres. The retail system for each centre may be 
considered to be open. 


The resulting boundary is defined in such a way that total 
expenditure available to the centres in the system must equal the 
total sales that result. 


In effect, a boundary is chosen such that the net boundary 
flow is negligible in comparison to the total trade in the study 
area. It is sited far enough from the centres under study to 
ensure that there is a negligible error in the calculations 
relating to these centres. 


The boundary is also drawn to include the majority of the 
people who are likely to do their regular shopping at the centres 
under study. 


There may still be a certain number of people living in areas 
outside the boundary who shop at the centres in question. An 
adjustment may be necessary to account for aberrations. If Total 
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SIMPLIFIED BLOCK DIAGRAM OF RETAIL GRAVITY MODEL 


Establish Consumer Establish Retail Estimate 
Preliminary 


Zones sa) Locations (3) 
Parameter 
Values 


(cee) 


Assembte Data on Assemble Data on 


Households, etc. Sar sh a0 OT, 


(C.) Shops, etc. 


J 


Estimate Distances or Travel 
Time from Each Zone to Each 
Centre from Grid Coordinates 


or Network (d; 5) 


Calculate Probabilities of 
Consumers in First Zone Being 
Attracted to Each Centre (P. .) 


VJ 


Repeat for Each Zone 


Sum Zonal Calculations for 
Each Centre 


Compare Calculated Sales with 
Actual or Estimated Sales for 
Each Centre (S.;) 


Adjust Parameter Values and 
Calibrate Model to Bring 
Calculated Centre Sales Into 
Accord with Actual Centre 

Sares Gos 6). A) 


Introduce Data on Introduce Data on 
Projected Future Proposed Expansion 
Populations and of Existingetentre 
Incomes and New Centres 


Introduce Data on 
New Centre Location 


Recalculate Market 
Shares Captured by 


Recalculate 
Market Shares 
Captured by Each 
Centre 


Each Centre 
Including New 
Centre 


Recalculate Market 
Shares 


Estimate Best 
Estimate Future Location, Size and 


Examine Impact of 
New Centres on- 
Existing Centres 


Retail Floor Planning of New 
Space Requirements Centre Development 


Fig. 2-3: Simplified Block Diagram of Retail Gravity Model 


McCabe, Robert W. Shopping Centre Dectstons: Evaluation 
Guides, Toronto, Province of Ontario, Ministry of Treasury, 
Economics and Intergovernmental Affairs, 1971. 


DETERMINING THE SPATIAL LIMITS OF A STUDY AREA 


1. The "Closed" System 


a) Flows to or from external areas subtracted from mass variables. 
b) External flows aggregated as a "dummy" variable. 


si 
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and E is external attractions 
e is external distance. 


2. The "Open" System 


a) Flows to or from external areas are aggregated to larger zones. 
b) External flows are calibrated through a series of rings. 


3 
Further interactions across 1. Iterate for all internal areas. 
the external boundary should PwecerateeroOr trips. tosdareas In 


be compensated for. LA era cee 
3. Iterate for subsequent trips to 
areas in 2, 3 and elsewhere if 
necessary. 


Fig. 2-4: Determining the Spatial Limits of a Study Area. 


Davies, R. L. and B. J. Styles (eds.). Gravity Models in Town 
Planning, Coventry, Lanchester College of Technology, 1969, 
p. 45. 


Available Retail Expenditures do not equal Total Actual Sales for 
the centres under study, then: 


(a) the study boundary is unrealistic and includes more people 
than those who do their regular shopping at the centres under 


study, or 


(b) a nearby large centre outside the study boundary is capturing 
a significant proportion of the retail sales from centres 
within the study area, or 


(c) there is a deficiency of shopping facilities in the study area 
and therefore a potential exists for a new centre or additions 
to established centres. 


Arriving at a study area boundary is less, likely to be a 
problem for a regional system of centres. The peripheral area is 
often devoid of competing centres. Problems do arise if an attempt 
is made to assess the potential of a number of centres that 
constitute only a portion of a metropolitan area. 


In this case a closed system is preferable and estimates may 
be needed to delineate the effective distance people might be 
willing to travel to patronize the centres under study -- the centre 
study boundary. 


Some account must then be taken of the centres existing outside 
the boundary to which a portion of the consumer spending within the 
centre study boundary might be attracted, especially from those 
consumers near the periphery. 


A dummy variable may be introduced to represent external 
centres and an appropriate distance measure assigned to it. (Fig. 
2-4) 


Consumer Origins 


The applications of the two retail models used in this study 
assume that the origin of all consumers is from their place of 
residence, that is, residential zones. 


The models have been capable of accounting for shopping trips 
made from place of work since their inception. Thus major 
employment areas can be properly accounted for as sources of sales 
for nearby shopping centres provided suitable survey data on 
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expenditures at place of work is available. 


The percentage of sales that come from the large number of 
office employees working in the central-city office buildings of 
major metropolitan areas like Toronto, Ottawa, and London provides 
a clear illustration of how significant this factor can be in some 
cases. Large decentralized industrial or office parks can have a 
Similar effect on suburban shopping centre sales. 


Adjustments in the average per capita retail expenditures in 
residential zones become necessary if place-of-work shopping is 
introduced. 


Centre Attractiveness 


Previous studies have shown that a variety of measures may be 
used as proxies for attractiveness. The number of items carried 
by the total complex was discussed earlier as a possible measure. 


Other measures have been devised which have proved to be 
equally reliable and less laborious in their preparation. Floor 
Space, annual sales, number of store types and number of employees 
are some of the alternative measures that have been applied. 


Shopping centre sales volume is probably the most precise 
measure. Unfortunately, information on sales by shopping location 
is not always available because of disclosure rules. For this 
reason the selling area in square feet -- gross leasable area -- is 
used either directly as a proxy or as the basis for estimating 
annual sales. 


Either option appears to be reasonably realistic for developing 
public policy because of the general consistency in sales and sales 
per square foot of gross leasable space in shopping centres. 


Total annual sales estimates from retail footage appeared in 
most trade area analysis examined in 1972 to be based on a 
relatively constant ratio of about $60 per square foot for 
Department Store Type Merchandise (DSTM) floorspace and $120 to 
$160 for Food floorspace. 


Separate estimates for DSTM and Food floorspace are usually 
preferable because of the significant difference in sales-per- 
Square-foot ratios. 


The limitations of these assumptions are fully recognized 
since it is known that sales-per-square-foot ratios can vary between 
centres and between the same stores in different centres. However, 
this is a problem that retail models share equally with trade area 


analysis. 


As a consequence, retail models cannot be held inferior in 
this respect. In fact, the need to be explicit in the assumptions 
made about sales at each centre indicate that retail models are 
more open than trade area analysis with fewer hidden assumptions. 


Floorspace figures expressed in square feet of gross leasable 
area are the most readily available measurement of attractiveness. 
There are no official published statistics on retail floorspace 
but this information can likely be obtained from planning offices. 


In many instances it is necessary to do an independent study 
of existing floorspace. Great care must be exercised in the 
methods of measurement and recording to ensure that the results 
are accurate and categorized in such a way that they relate 
exactly to published sales-per-square-foot data. 


Distance Measures 


Distance may be measured in several ways. The simplest is 
straight-line distance. The most complex is probably an index of 
time/distance/cost for each of the available modes of transportation. 
Other simple measures may be expressed in terms of travel time or 
travel cost. 


Straight-line distance takes no account of such matters as 
traffic congestion, the frequency and price of public transport, 
variation in car ownership, cost of parking, physical barriers and 
standard of roads. It seems clear that some combination of time- 
distance and travel cost provides a preferable measure. 


In certain situations, however, very little in accuracy 
appears to be sacrificed in using straight-line distances. In 
other cases it seems to result in significant distortions. 


Straight-line distance has the advantage that it can be 
computed internally by the computer using a grid reference system 
saving considerable time and expense. However, the use of straight- 
line distances in large metropolitan areas can usually be avoided 
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if good travel-time information at a fine scale of detail is 
available and if the time and cost involved in retrieval from 
traffic study data files is not exorbitant. 


In smaller centres, where traffic data is less likely to be 
available, the use of straight-line distance presents more serious 
limitations because small centres often have a less homogeneous 
transportation system and a more diverse settlement pattern. In 
some circumstances it is possible to actually drive the distances 
to all centres from each population zone in the study area under 
varying traffic conditions and thus produce a valid set of travel- 
distance figures. 


The two major assumptions about the effects of distance on 
shopping movements are that they are the same for all socio-economic 
groups within the study area and that the parameter value of the 
distance variable remains constant over limited time periods. 


Future populations are viewed as reacting in the same way to 
the deterrence factor of travelling as in the base year, notwith- 
standing possible effects of changes in, say, disposable income, 
increased leisure, improved transportation facilities or altered 
marketing strategies. The effect of this assumption must be 
weighed against the observed dramatic change in consumer food 
shopping patterns that occurred in 1970 when the leading chain food 
Supermarket went to discount pricing in Ontario (Tigert 1971). 


New travel-time estimates may be used to reflect road or highway 
improvements but parameter values cannot be adjusted with any degree 
Of assurance. Projections limited to five years into the future 
may partially offset this problem but short-term changes in 
consumer behaviour cannot be taken into account. 


New Centre Locations 

Retail models do not automatically select optimum location for 
new centres. Alternative sites or arrangements of sites may be 
evaluated against whatever criteria have been judged as desirable. 


Possible shopping locations are first decided upon. Estimates 
Of probable sales levels at the new centre and all the existing 
centres are then generated by the model. The consequences of other 
contemplated actions can be estimated by analyzing alternative 
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locations and by varying the size of proposed centres. The 
solution finally adopted will be optimal only for the range of 


alternatives tested. 


Several types of analysis are possible. The analyst can 
examine the minimum size needed for a proposed centre, the possibility 
of adding more retail floorspace in any location up to a specified 
maximum, or the effect of alternative population densities around 
the centre. The analyst's imagination and the realities of the 
retail situation are the only limits on the analyses possible. 


Both models offer considerable flexibility to the planner. The 
Huff version allows for the use of consumer survey data to calibrate 
the model where such data is available at a zonal level. It also 
allows for the evaluation of the profit potential at a site under 
incremental size levels up to a specified maximum. The Cullen 
version also has a number of options that can be selected -- a 
single-run projection, an equilibrium projection and a dynamic 
projection. 


All procedures for both versions were tested during the course 
of the study. While no attempt is made here to discuss all the 
options and their applications, it was confirmed that they do 
operate as claimed. Nonetheless, further research in the use of 
the options in operational applications in Ontario would be helpful, 
and instructions for implementing the various options are included 
in the Technical Appendix. Further discussion of the options will 


be found in the Huff and Cullen references. 


Retail modelling has qualities that appear to make it superior 
to other methods. There are, however, a few factors affecting 
locational choice that are not considered. These include site 
characteristics such as visibility, total road frontage, condition 
of road access, topography, drainage, amount of grading required, 
soil conditions, services available, cost of purchase, municipal 
building restrictions, and performance standards. 


These factors are important to planners and retailers. But 
they are not critical since it can be assumed that the sites to be 
evaluated would be expected to meet these conditions in any case. 


SOME GENERAL COMMENTS ON RETAIL MODELS 


The principal criticism of retail models has centred around the 
impression that the retail gravity model is based on a weak analogy 
between aggregate consumer behaviour and Newtonian physics. 


Early versions, like Reilly's Law of Retail Gravitation, were 
indeed dependent on the similarities between the attractiveness of 
mass in retailing and the deterrence of distance. 


New strength and new direction has been given to retail modelling 
in the past two or three years by major conceptual and Operational 
breakthroughs. 


This is particularly true of the work of Wilson at the Centre 
for Environmental Studies in London and that of the scholars and 
practitioners working at the Planning Research Applications Group 
in Reading (Wilson 1969c and Wade 1973). 


The Wilson derivation, while resembling in general form the 
older gravity models, overcomes any possible criticisms of the 
retail model as having weak empirical or conceptual bases. The 
Wilsonian entropic model used in these studies was rigorously 
derived through the process of statistical mechanics. It is based 
on the principles of entropy maximization which hypothesizes that 
the most likely state of any system can be approximated given 
sufficient knowledge about how the system behaves. 


If all that is known about the system is reduced to a series of 
System equations, the most probable state of the system can be 
calculated. In the case of the macro-retail shopping system this 
has been accomplished and remarkable replication of observed 
Shopping activities results. 


An important feature of the Wilsonian version of the model is 
the use of an exponential form of the attractiveness and deterrence 
functions in place of the power functions employed by Huff and 
Lakshmanan. 


The shape of the exponential curve of parameter values is 
Somewhat less flat than that of the power curve. The calibration 
of actual to calculated results is therefore more easily obtained 
as the differences between solutions are more sharply defined along 
the curve of exponential values. 


The effect of the exponential function is to decrease the number 
of longer trips hypothesized by the Huff version in favour of shorter 
trips, thus producing a distribution more nearly like that typically 
observed. 


A criticism of the Cullen version is that, like the Lakshmanan 
version, it calibrates on two-parameter Values. stl hemtechnicalieaies 
of the propriety of the two-parameter version are the subject of 
much discussion in the journals (Batty and Saether UG 720s 
Lakshmanan overcame the difficulty by setting the attractiveness 
parameter at 1.0 and calibrating only on the distance parameter. 
This simply produced a restatement of the Huff version. 


A more serious problem is that of determining whether the 
calibrated parameter values provide a unique solution or whether 
other values may not provide an equally good fit of actual and 
calculated values of centre sales. 


This problem has not yet been resolved. However, if due care 
is taken to ensure that the area of search for a good fit excludes 
the "bogus" calibration that is achieved when both values are set 
equal to zero, and that it includes values well beyond those 
typically found in studies conducted elsewhere, a set of values 
will emerge that appear to be optimal. 


In this study, and in others conducted before and after, a 
best fit was found within narrow limits in each case that optimally 
related actual and calculated centre sales. 


The most frequently asked question about retail models is why 
they fail to take into account all the factors that are assumed to 
affect centre sales. In fact, versions of the model have been 
formulated in which other measures of attractiveness in addition to 
size have been incorporated. 


For example, it is known that the reputation of a particular 
retailer may have a beneficial effect on the sales in a centre vy 
that particular retailer located in the centre. That retailer's 
drawing power relative to size is larger than that of other 
retailers. 


The benefits of the selected retailer's presence presumably 
will be reflected in the operational efficiency of the centre and 
consequently in the level of sales per square foot derived from the 
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selling space. If these beneficial effects could be measured 
without violating the disclosure rules of Statistics Canada, great 
advances would be possible in retail modelling. 


In these circumstances, given the mix of shops in a centre, the 
attractiveness parameter might be estimated independently of the 
model and a valid calibration achieved on the distance parameter 
alone. Other differentials in drawing power by class of retailer, 
by geographic areas, and by the socio-economic characteristics of 
customers can be envisaged. 


This exact information is contained in the files of Statistics 
Canada and could be published subject to the normal disclosure rules. 
A shortage of government funds apparently prevents this being done. 


Most other questions can be answered by stating that retail 
models use the same data sources as trade area analysis for sales, 
centre sizes, populations, incomes and retail expenditures for the 
past, present and future. 


Both can encompass the results of consumer survey if necessary. 


Shopping models, however, incorporate spatial relationships in 
a much more precise and rigorous manner and allow for the testing of 
a large number of alternatives which trade area analysis cannot do, 
except at prohibitive levels of expenditure of time and money. 


Retail models are macro-level in nature. They do not attempt 
to predict the behaviour of individual consumers although there is 
little doubt that the understanding of observed shopping patterns 
would be enhanced by a clearer understanding of the shopping 
behaviour of individuals. 


Such micro-level models are now being studied and much has been 
learned. Still, it is questionable, even if it could be explained 
today why two families of identical socio-economic and cultural 
backgrounds living side by side consistently shop in different 
shopping centres, whether any better prediction of consumer shopping 
behaviour at a macro-level could be achieved. 


The real issue for urban and regional planners is not whether 
Shopping models might be improved, but the degree to which planners 
themselves are attempting to apply what is already known when 
answering the complex questions with which they are faced. 
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At the macro level alone, a considerable amount of 
preparation is required before shopping models can be used in the 
decision-making process. One essential requirement is a body of 
reliable statistics, at a fine level of detail, on such fundamental 
matters as retail and service floorspace, population distribution, 
distance and income estimates, retail spending and sales. 


Just as essential is a sound grasp of the nature of retailing 


as an urban activity. 
CONCLUSIONS 


Whatever their real shortcomings, retail gravity models have 
been developed in the past decade into a powerful technique for 
interpreting the behaviour of aggregations of people at a macro 


scale. 


They still appear to be the best means available to date for 
evaluating a number of shopping centre proposals at one time. 


They represent the only methodology aimed at determining the 
impact of new centres on the sales of all other existing shopping 


facilities in a municipality. 


The flexibility of the retail model for testing a wide variety 
of alternative patterns of retail location remains largely 
unchallenged. 


The gravity model's open-system approach and its underlying 
principle of choice behaviour makes it unnecessary to arbitrarily 
derive an exclusive trade area for each individual centre. 


Of the methods currently developed for studying retail area 
structure, the concepts of the retail model are the most sensitive 
to the problems of bridging the gap that exists between strictly 
empirical analysis and the development of a true model of individual 
consumer choice behaviour from which parameter values may be 
estimated directly. 


Huff's introduction of a probabilistic consumer-choice- 
behaviour formulation was a significant step forward. It approached 
the problem of evaluating and predicting the size of retail centres 
in a given locational pattern from a consideration of the consumer's 
spatial behaviour. 


Huff acknowledged that consumers do not appear to inexorably 
patronize the nearest centre by developing a conceptual framework 
in which it is possible to assess the probabilities of choice made 
by consumers among a series of competing centres varying in size 
and separated by unequal distances. 


The notion of non-uniqueness in consumer shopping patterns 
differed significantly from all previous approaches. Both the 
subjective evaluation methods of Central Place Theory and trade 
area analysis, and the early retail gravity models explicitly or 
implicitly assumed exclusive trading areas. Instead, Huff 
described trade areas in terms of a series of radial zonal 
probability contours, and not as a fixed line circumscribing a 
shopping centre. 


The probability formulations of the retail model offer a 
significant improvement over earlier models. They replicate the 
interaction between a continuous population distribution and a 
large number of shopping centres representing change from a two- 
centre interaction model to a many-centred interaction model. As 
a consequence, there has been a change in emphasis from the market 
area of the seller to the shopping area of the buyer. 


The Huff, Lakshmanan and Hansen, and Cullen versions all permit 
the evaluation of a system of commercial shopping centres at one time. 
They make it possible to conduct policy experiments by adjusting the 
variables within the model to discover instantaneously the effects 
of various alternatives and assumptions. These features make the 
method a very valuable tool for planners and an important source 
of information for government decision makers. 


A great deal of exciting new work has taken place in England 
on extensions of the retail model to incorporate all levels of the 
Shopping hierarchy and the multiple groups of trading areas found 
in large poly-centred urban regions. Much of this material has now 
been made available (Wade 1973). 


Further progress in testing these novel applications in Ontario 
Can be expected. The implications for resolving retail location 
problems at a municipal and regional level are readily apparent, 
particularly for those concerned with developing planning policy 
guidelines for the orderly growth and development of retail 
facilities within newly created regional governments. 
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Further details of the options which are available to the 
planner are discussed in the case studies chapter. Instructions on 
their use can be found in the Technical Appendiz. 
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3 The Case Study Municipalities 


Purpose and Selection of Municipalities 


Applications of retail modelling were made in two Ontario 
municipalities during the course of the study. Both the Huff and 
Cullen versions were tested in each municipality. 


The purpose of the applications was to see to what degree the 
models were able to provide answers to the sets of planning and 
policy questions posed in Chaprer. i 


The case study municipalities were selected on the basis of 
their size, complexity, availability of information, and the degree 
of interest shown by the staffs of the planning boards concerned. 


Belleville was one of the centres chosen for study. 
Considerable information about the retail trade area of Belleville 
was already available from Shopping Centre Decisions: Evaluation 
Gutdes, and from a trade area analysis conducted in 1969 prior to 
the opening of the Quinte Mall. Thus a convenient basis was 
available for testing the predictive Capability of the retail 
models. 


The second location selected was the area of the Central 
Ontario Joint Planning Board. While not a municipality in the 
Strict sense, it was apparent at the time of the study (summer 1972) 
that the area under the Board's jurisdiction would be central to 
whatever boundaries were eventually drawn for the Regional 
Municipality of Durham announced at the time and now established. 


The Board had recently completed a full inventory of retail 
and service floorspace. Moreover, the Oshawa Area Planning and 
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Development Study provided a rich source of data on population and 
economic growth potential. The planning board staffs of both 
municipalities were keenly interested in participating in the study. 


Study Approach 


A general hypothesis was formulated about the capabilities of 


the retail model: 


Retail models are capable of evaluating a wide 
range of alternative arrangements for the 
provision of retail shopping facilities under 
varying assumptions about future travel times 
and distances, levels of retail expenditure, 
and population growth rates, distributions and 
densities. 


The validity of this postulation was tested through a number of 
working hypotheses related to the specific capabilities of the models. 
It was possible to make comparisons between models and between case 
studies by testing the Huff and Cullen versions of the retail model 
in the two case study municipalities over a number of time periods. 


Here are the seven working hypotheses that formed the basis 
for the testing conducted in each municipality: 


ie Measures of attractiveness, distance, consumer expenditure and 
shopping centre sales are available for Ontario municipalities 
in enough detail and at an acceptable standard of reliability 
to allow the models to operate. 


(de The models, when calibrated, provide a reasonable replication 
of.existing patterns of retail trade in the case study 
municipalities. 


or The calibrated models have a short-term predictive capability 
sufficient to make valid projections for at least five years 
into the future and for longer periods where an adequate basis 
is available for predicting all the factors that enter into the 
models. 


4. The models indicate the probable consequences of introducing 
an added centre or an expanded centre into the municipality. 


oe The models produce a useful comparative measure of the levels 
of service to customers provided by alternative possible 
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locations. 


Gs The models produce a useful comparative measure of the 
productivity for retailers provided by alternative locations 
and alternative sizes of centres. 


7h The levels of customer flow and expenditure flow between 
residential areas and shopping centres closely approximate 
known conditions. 


The Use of Terms 


It is often difficult in studies of this nature to clarify the 
distinctions between planning boards and their staffs, and among 
the legally defined boundaries of the planning area, the 
municipalities concerned, and the trade area that has been defined 
for study purposes. 


In this chapter, the municipal agencies with which the study 
team dealt were the staffs of the relevant planning boards. For 
simplicity in presentation the distinction between staffs and 
boards has been abandoned and only the board referred to. In 
Oshawa, for example, it was the staff of the Central Ontario Joint 
Planning Board who were consulted. The report refers simply to 
Oshawa or the Board. Area is used as a reference to the planning 
area of the COJPB as adjusted for purposes of the study. 


Similarly, Belleville is used when referring to the trade area 
boundaries surrounding Belleville as defined in the study and not 
to the boundaries of the Belleville Planning Board. In actual 
fact, the study team worked with the staff of the Quinte Planning 
Board. 


Moreover, the boundaries of the Quinte Planning Board's 
planning area do not correspond to the trade area boundaries of 
Belleville as used in this study. Accompanying maps have been 
provided in the text to make sure that the boundaries used for 
purposes of defining areas for this study are clear (Maps 3-2, 
3-3 and 3-4). 
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THE OSHAWA CASE STUDY 


Sources of Data 


i Study Area 

The area of market influence of the retail facilities in the 
Oshawa area corresponds closely to the boundaries established for 
the Central Ontario Joint Planning Board. The boundaries 
established for the study were accepted as those of the Board. 
They encompass a relatively self-contained regional trade area 
influenced by Metropolitan Toronto to the west (Fig. 3-1). 


rae Zones 


The Oshawa area had been divided into 67 Basic Planning Units 
of varying size. The smaller sized units occur in the more densely 
populated urban areas and the larger sized units in the rural areas. 
The number, size, and population in each of these units justified 
their use as population zones for the study (Fig. 3-2). 


zs Centres 


An extensive commercial inventory was completed and published 
for the study area in the spring of 1972. As a consequence, 
shopping centre floorspace information was available at the 
required level of detail. 


A review of the survey established 45,000 square feet of 
leasable area as the smallest size for a centre to be included in 
the study. The eight centres over this size provided a range of 
comparison and shopping goods large enough to be influential at a | 
community or regional level. 


Where the floorspace was below 45,000 square feet, the 
predominant proportion of the total floorspace was devoted to food. 
Centres-in this lower category appeared to have a strong local | 
influence but did not seem to have a significant impact at a | 
regional scale. 


4, Floorspace Estimates 


Floorspace estimates for 1971, both for Food and DSTM 
(Department Store Type Merchandise), were taken from the Board's 
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1972 report on commercial structure. Personal service floorspace 
was deleted from both the planned and unplanned centres. It was 
found that the majority of personal service floorspace was located 
within the small unplanned centres not included in the centres to 
be studied. 


Floorspace figures for 1966, 1976 and 1981 were provided by 
the Board. The 1966 figures were prepared by examining the opening 
Or expansion dates of existing centres and recreating the 1966 
centre sizes and locations. 


The figures for 1976 and 1981 on possible future centres, and 
their number, location, and floorspace, were confined to those that 
had been proposed by the constituent planning boards in the area or 
discussed by private developers with them (Table 3-1). No other 
detailed examination of future retail floorspace requirements was 
undertaken. 


ie Zonal Population 


Zonal populations for the study years 1971 through 1981 were 
expressed in terms of households. As part of a planning information 
system, the Board had collected information by households for all 
67 zones for the year 1971. Estimates of the total number of 
households in 1976 and 1981 were based on planning studies and 
growth projections for the area. Distribution among zones was 
provided by the Board through reference to local official plans 
and other planning documents in their possession. 


oF Distance Estimates 


The Board prepared a travel-time matrix from recent traffic 
studies simulating the 1971 situation. Times were indicated for 
all junction-to-junction links covering routes from each zone to 
each centre. It was assumed that travel times from the zone 
centroids to centres would remain unchanged up to 1976. This 
assumption was later extended to 1981. 


The use of straight-line distances in place of a travel-time 
matrix is possible in the Huff version of the model. A test of 
this aspect of the model was considered to be valid in Oshawa 
because of the more uniform size of the population zones. Straight- 
line distance measures were computed from each zone to each centre 
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internally to the Huff program from grid co-ordinates of the 
centroids of each zone and each centre determined by the study 


team. 


The population centroids of zones are likely to shift over 
time. It would be preferable to take these shifts into account 
but it was not possible to do so in this study. 


73 Income Estimates 


The derivation of 1971 income information by household for each 
of the 67 study zones proved to be a very difficult procedure. Some 
income information was available from earlier transportation studies, 
from other population information known to the Board, from federal 
tax statistics, and from the knowledge of local and provincial 
agencies and informed officials. 


A random sample of family incomes for selected zones was also 
undertaken as a means of validating the estimates for 197] through 
the co-operation of the Area Credit Bureau. The final average 
household income figures for each zone were arrived at from the 
Survey by applying a family/household ratio. 


8; Food and DSTM Expenditures 


The proportion of household income available LOnere banda 
spending was derived by reference to Family Expenditures in Canada, 
1969. The percentage of total income spent on retail goods was 
calculated for each level of income. The total income per household 
for each zone was then adjusted by the approximate percentages to 
estimate Food and DSTM expenditures per household. 


The Financial Post's Survey of Markets 1971 reported that 
Personal Disposable Income for Oshawa was 16.3 percent greater 
than the Ontario average while retail spending was only 10.6 
percent above the provincial average. 


These Figures suggest that there is a five to Six percent net 
outflow of potential expenditures on retail goods. It was estimated 
that 100 percent of the potential Food expenditures are made within 
the study area. As a consequence, the income available for DSTM 
purchases in each zone had to be reduced enough to ensure that 
total expenditures would equal total sales. 
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The final Food and DSTM zonal figures for 1971 were then 
increased by 1.5 percent per annum to account for increases in 
1976 and 1981 expenditures in real dollars. A 5.7 percent net 
outflow was held constant over the study period. 


9. Estimates of Centre Sales 


Further reductions in available Food and DSTM expenditures per 
household were required because the group of centres under 
consideration represented less than the total available shopping 
floorspace. 


Food and DSTM floorspaces for the eight centres were expressed 
separately as ratios of the total Food and DSTM floorspace in the 
area. The estimated average per household expenditures on Food and 
DSTM were then reduced by these ratios for 1971, 1976, and 1981. 


Estimated actual sales for the eight centres in the area were 
calculated by multiplying the Food and DSTM floorspace in each 
centre by estimates of sales per square foot varying from $110 to 
$120 for food and from $50 to $60 for DSTM. At levels of sales of 
$120 and $60 respectively, the estimated actual sales for the eight 
centres were $93,500,000. 


When the households in each zone were multiplied by the 


corresponding estimated 1971 Food and DSTM expenditures per household 


the calculated sales potential was approximately $93,000,000. 


CONCLUSIONS - DATA SOURCES 


The data needs of the retail models are extensive. In Oshawa 


it was found that there was little difficulty in establishing a 


suitable set of trade area boundaries, in defining statistical 
areas suitable in number and size for use as population zones, in 


delineating centres, and in measuring the appropriate floorspaces, 
zonal populations and travel-time distances. 


Difficulties arose, however, in measuring the retail 


expenditures of households in each of the population zones and 
estimating the actual sales of shopping centres. 


Fortunately, the models are constrained in such a way that 


total sales in the centres included must equal total expenditures 
available. By carefully estimating each total separately, a 
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reasonable measure of the validity of the estimates can be obtained. 


In Oshawa, for example, the difference between total expenditures 
calculated from estimates for 67 zones to total retail sales derived 
from floorspace measures and sales per square foot was less than 


half of one percent. 


Part of the reason for the availability of usable data in 
Oshawa was the Board's awareness of the fine level of detailed 
planning information that is required to support regional planning 
policies and to measure their impact on small zones within the area. 
Problems of estimating zonal incomes and expenditures are found in 


many planning studies. 


Similar problems of approximating actual sales for individual 
centres do not appear to have been resolved elsewhere. Both kinds 
Of problems are common to trade area analysis and retail models 
although retail models appear to have the advantage of requiring 
that detailed checks be made on the validity of the estimates 
entering into the program. 


It can be concluded that Oshawa had adequate information 
Sources available to provide the data required to operate the model. 
Additional research on means of estimating incomes and expenditures 
for small statistical areas and on better ways of approximating 
actual centre sales is indicated, although the methods used in 
Oshawa appear to be more refined than those found in comparable 
studies. 


Options Tested 


The two versions of the retail model were tested using Oshawa 
data. The decisions on which sets of options to test were made in 
close collaboration with the Board. The Board subsequently asked 
the study team to incorporate only selected options. 


Their suggestions included considering which existing retail 
centres might expand, which of the centres might be operational 
by 1976 and 1981, and which combinations of centres and sizes 
seemed most reasonable. 


The option decisions were, in effect, an attempt by the board 
to anticipate or simulate the future retail structure within the 
Study area from information available to them. 
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the "bogus calibration" characteristic of the process. When X = 


’ 


Table 3-1 provides a summary of all the options tested and the 
corresponding total floorspace. 


Map 3-2 depicts the location of each retail centre listed in 
Hable 3-1. It ts worth noting the relationship of the centres to 
major transportation corridors through the area and, thus, to major 
shopping complexes on the eastern periphery of Metropolitan Toronto. 


The results of only two options beyond 1971 are discussed here 
although the results of tests of all options were made available 
to the Board. 


A total of twelve shopping locations were considered in 1976 
- Option 3 of which Valiant Property Management (Oshawa), Western 
Auto Parts (Oshawa), F & T Developments and the S. W. Corner of 
Dundas West and Francis St. (Whitby) are additions to the 1971 
inventory. 


The sizes of the original eight centres were held constant for 
1976 - Option 3, no growth or reduction in retail floorspace for 
those eight centres. 


The 424,700 square feet of retail space added by 1976 represents 
a further decentralization of retail activity within the area. But 
the largest share, 63 percent of the total added, was situated in 
Oshawa itself. 


A review of 1981 - Option 3 suggests the POssal. Duby Ot. a 
further increase in the amount and dispersion of retail floorspace 
within the study area. Fourteen retail centres are included 
representing an addition of two new centres over the 1976 level. 
The two new centres, Riznek Construction (Darlington Township) and 
Bradley Farm (Whitby) provide for a retail floorspace increment of 
307,000 square feet. The floorspace levels of the other twelve 
centres were held constant beyond those shown in 1976 - OptiLOn: 3. 


Calibration 


The Huff version was calibrated on 1971 data by varying the 


“parameter value from 0.0 to 2.0. With A = 0.0, the model produced 
a nearly perfect fit between calculated and actual sales confirming 


OF75, an acceptable fit was obtained. The fluctuations between 


calculated and actual sales on either side of the calibrated value 
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were less pronounced than those observed in the Cullen model. The 
centres were appropriately ranked in order of their sales magnitudes. 


The Cullen version was calibrated on 1971 data varying the 
value of @ in:steps of 2001) ‘and ‘the X ‘value “Injstepsiof 7005 aaa 
minimum in the summed difference squared between calculated and 
actual sales occurred when a was set at 0.995 and A at 0.011. The 
values on either side of this calibration point produced fits 
between calculated and actual sales that were markedly different. 


The summary site information for Oshawa for 1971 at the 
calibrated value of lambda (A) using the Huff version is shown in 
Table 3-3. The comparable summary for the Cullen version is shown 
in Table 3-4. Differences can be accounted for by the fact that in 
Table 3-3 an average of $72.00 per square foot was applied to all 
centres to estimate the actual sales. In the Cullen version 
Separate sales-per-square-foot values were applied to the Food and 
DSTM components to estimate actual sales. A review of calculated 
and actual sales suggests that the Cullen version gave a better 
replication of what is estimated to be the pattern of retail 
shopping in the Oshawa area. 


Table 3-5 is included to show the effects on summary site 
output when straight-line distances have been calculated internally 
using grid co-ordinates in place of a travel-time matrix. A 
comparison of Table 3-5 with Table 3-3 reveals only slight 
variation between the calculated and actual sales levels for 
individual centres. 


It might be suggested that in the Oshawa situation there is 
little difference regardless of which measure is used. However, 
the significance of the differences has not been established. It 
can be assumed that if straight-line distances were used throughout 
the calibration, some small adjustment in parameter values would 
have resulted. 


It can only be concluded that the use of either straight-line 
or travel-time distances should be consistent within a given study. 
Further research on this problem is indicated. 


Predictive Capability 
The predictive capability of the model was tested by holding 
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the calibrated parameter values constant and attempting to 
reproduce a known situation five years ahead, a forecast, or five 
years behind, a backcast. In Oshawa, a backcast was undertaken 
from 1971 to 1966 -- a prediction of the 1966 situation using 1971] 
parameters and 1966 data. 


Tables 3-6 and 3-7 illustrate the results of the backcast for 
the Huff and Cullen models. The Huff backcast was adequately 
reliable for the largest centre, but relationships and rankings 
were lost for the remainder. In the Cullen backcast, when allowance 
is made for the difference between total synthesized sales and 
actual sales, the goodness of fit is remarkable. 


It can be concluded that for Oshawa the Cullen version confirmed 
the short-term predictive Capability of the retail model. The Huff 
version was somewhat less reliable. It was apparent that the Huff 
version would be improved if the program was adjusted to allow 
separate sales-per-square-foot estimates to be applied to the Food 
and DSTM floorspace components of each centre when estimating actual 
Sales. This adjustment has now been made. 


Impact of New Centres 


Tables 3-8 and 3-9 show the results of the use of the Huff and 
Cullen versions to evaluate proposals for additional centres to 
meet the need for added retail floorspace arising from the increased 
household and the increases in the level of real income per household 
meayected in this study. 


The arrangements proposed by the board to meet the needs created 
by a sales increase of about $54 milion. trom $93 million in 1971. 
(Tables 3-3, 3-4) to $147 million in 1976 (Tables 3-8, 3-9) fall 
Short by about $22 million -- the difference between calculated and 
actual sales. 


At a level of $72 per square foot average for Food and DSTM 
space this would suggest that there is room for an additional 
300,000 square feet of retail Space by 1976 beyond that suggested 
Miet276 - Option 3 (Table 3-1). 


In effect, two million square feet of added retail floorspace 


will be required between 1971 and 1976 if the predicted population 


increases and income levels are matched, and if consumers are to 
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continue to have comparable levels of retail shopping facilities 


available to them. 


The underestimation of the amount of retail floorspace required 
to meet future needs could have been overcome if the model had been 
used as a means of projecting total space needs. This requires 
special reliance on the household and income projections. In Oshawa 
these appeared to be sufficiently robust that this use of the model 
might have been explored if resources had permitted. 


Table 3-9 also indicates that the largest decline in market 
Share will occur in the Oshawa Centre, from 35.49 percent to 26.27 
percent. It would be interesting to rerun both models to see what 
the consequences would be of adding the required additional 300,000 
Square feet in the City of Oshawa or elsewhere. 


The fact that the option table didinolLeprovidessursicient 
additional floorspace is reflected in the sales-per-square-foot 
figures. Calculated sales per square foot climb to about $13 above 
the level that would be expected if retail productivity per square 
foot remained at the 1971 levels in real dollars. 


Similar forecasts were made for 1981. 1981 - Option 3 indicates 
a shortfall of $70 million in real dollars over the 198] requirements, 
reflecting the need for one million square feet of retail space 
beyond what is provided for in that option (Tables 3-10, 3-11). 


An analysis of Tables 3-10 and 3-11 shows a continued drop in 
the market share for established centres as new centres are opened, 
and a continued rise in calculated sales per square foot as a 
result of major under-storing of the area. 


Levels of Service 


Two measures of level of service are available. The first -- 
the level of sales per square foot generated -- is provided by both 
versions of the model. It can be assumed that if the sales per 
Square foot of calculated sales for a centre are beyond what could 
be reasonably expected, say, $110 per square foot compared to an 
average for the area of $72, the consumers in the area of the 
centre are not being well served. Conversely, a less than average 
sales-per-square-foot figure might indicate an overstored condition 
in a local area. 


The second measure is provided only by the Cullen version in 
the column marked AVTL, or average trip length. When comparing two 
alternative arrangements of retail floorspace to serve an area 
population, one measure of service to customers is the average 
length of the shopping trip for all customers. 


In Table 3-4, for example, the average trip length for the 
shopping pattern represented by the synthesized sales is 10.08 
minutes, ranging from a low of 8.32 minutes for Centre 4 to 15.49 
minutes for Centre 5. If alternative arrangements of centres are 
M@esible -- say, the three options in 1976 or 1981 in Table 3-1 -- 
the differences in average trip length among the three options could 
be one basis for selecting one option over another. 


Comparative Measures of Retail Productivity 


The Oshawa study did not provide a suitable means of 
illustrating the usefulness of comparative measures of retail 
productivity because the options suggested for testing did not serve 
the area adequately. Other uses of sales per square foot as a 
measure have been mentioned, however, in two previous sections. 


Customer and Expenditure Flows 


Each of the two versions of the retail model produces 
information on the amount of retail expenditure available in each 
zone and on the flows of customers and expenditures from each zone 
to each centre. 


In the case of the Huff version, this information is provided 
in the form of neighbourhood summaries (Table 3-12). The information 
that has been provided to the model about the zone is summarized in 
the first line. A table shows the number of customers, the amount 
expended in each centre, and the total zonal expenditure as well 
as other information. 


ine the case of the Cullen version, similar information 1s 
provided in a series of tables. Zonal expenditures are provided as 
a single table (Table 3-13). Zone-to-centre customer and expenditure 


BOWS are shown aS matrices (Tables 3-14, 3-15). 


The matrix version provided by the Cullen version is more 
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compact and easier to read. The Huff version provides a neat set of 
tables but the sheer bulk of the printed results and the added 


costs in computer time limit their use. 


If consumer shopping survey data were available it would be 
possible to produce two comparable matrices of customer or 
expenditure flows -- one from survey and one from the output of 


the calibrated model. 


The possibility of doing so inethe counsesorethissrescarchamac 
discussed but the costs of the survey were beyond the resources 
available to the study team. Comparison of the results would do 
much to confirm the reliability of the retail model. Hopefully, 
this will be done in future applications. 


Conclusions - Oshawa 


The conclusions on data needs were discussed earlier. The 
options suggested for testing did not provide enough additional 
retail floorspace to match the forecasted needs. This made analysis 
of the results of testing less rewarding than anticipated. 


Both models calibrated easily on the 1971 base year data. The 


backcast to 1966 confirmed the predictive capability of the models, 
with greater reliability being shown by the Cullen version. 


The shortfall in the total amount of floorspace included in the 
option to be tested confirmed the potential value of the models for 
estimating future floorspace needs. 


On the basis of the Oshawa study, the Board was provided with a 
quick analysis of their initial#@intuitive projections and the means 
to progressively refine their proposals to a smaller number of 
feasible solutions. The potential usefulness was shown of measures 
of service to households, retail productivity, and customer and 
expenditure flows. 
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Belleville were based on the Board's knowledge of local conditions 
and on a detailed market analysis done in 1969 by a private 
consulting firm (Map 3-3). By choosing boundaries identical to 
those used in the 1969 study of Belleville's market potential, it 
was possible to compare the results of the trade area analysis and 
the retail model. 


The adopted boundaries can probably be considered conservative 
by the time that the 1981 projection year is reached. The drawing 
power of the estimated facilities needed to serve the study area 
population in 1981 is likely to extend beyond the defined boundaries, 
particularly to the west and north. 


a. Zones 


Forty-four wards and subdivisions of wards were adopted as 
population zones within the corporate boundaries of Belleville (Map 
3-4). Ward subdivisions constitute the city's polling districts. 
They are the smallest unit for which practical retrieval of 
population data could be achieved by the Board. 


Outside the City of Belleville 13 very large zones were 
established to cover the remainder of the area. These were made up 
of single urban centres or a combination of townships and smaller 
municipalities comprising the three large area units defined in the 
1969 study, Secondary Trade Area, Tertiary Trade Area A, and 
Tertiary Trade Area B. The concentration of so many zones within 
Belleville is questionable. A more equal distribution of the number 
of zones between the inner and outer areas would have been preferable. 


S.. Centres 


A detailed inventory of existing commercial retail facilities 
Was available from the 1969 study. The Board was able to arrive at 
data for 1966 and 1971 by a careful review of centres opened before 
and after the 1969 inventory. The majority of major centres in 1966 
and 1971 were unplanned. A minimum floorspace of 50,000 square feet 
Was adopted to ensure the general mix of store types and goods in 
the centres indicated that a community or regional function was 
being performed. 


4. Floorspace Estimates 

Floorspace estimates for both Food and DSTM were developed from 
the 1969 inventory. For the unplanned centres those groups of Food 
and DSTM stores were included which fulfilled the minimum criteria 
defined in the consultant's report. Only stores that did meet the 
standards were felt to contribute to the attractiveness of the 


unplanned centres. 


The size and number of possible future planned centres and 
their location were obtained through discussions with the board 
(Table 3-2). Some centres included for consideration arose from 
past proposals. Some were proposals that were in the early stages 
of discussion with private developers. Others were of interest to 
the Board as a means of evaluating alternative patterns of growth 


and development. 


The Food and DSTM floorspace of potential future centres was 
derived from preliminary submissions made by private developers or 


from estimates made by the Board. 


5. Zonal Population 


Household information for Belleville was not readily available. 
In consequence, 1966 and 1971 populations for the 44 zones within 
Belleville were derived from Statistics Canada reports and local 
data files. The projections for 1976 and 1981, and their 
distribution within the corporate area, were supplied by the area 
planning board using official plan references and a labour force 
index-population ratio. 


Populations for 1966 and 1971 in the 13 outer zones were taken 
from the provincial municipal directories. 1976 and 1981 population 
projections were based on estimated ratios reported during the 
1969 market study. Independent estimates were rejected in order to 
maintain as great a degree of comparability as possible between the 
trade area study and the retail model results. 


6. Distance Estimates 


A travel-time matrix was not available for Belleville. As an 
alternative, straight-line distances along major roads from each 
zonal centroid to each centre were measured in decimal inches from 
two differently scaled maps, one for Belleville and the other for 
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the peripheral areas. 


This required a scale conversion procedure. Converted figures 
were then added to indicate the full distances between the zones and 
centres on the two map schedules. The procedure here accommodated 
the many small zones within Belleville and the few macro-scale 
zones in the very much larger outer area. The resulting measures 
were employed as actual distances in the tests of both 
versions. 


Although the Huff version of the model has the capability to 
calculate straight-line distances internally, given grid co-ordinates 
for zone centroids and centres, the differences in scale of the two 
maps used effectively precluded this option. 


To help overcome the limitations imposed by the use of 
actual distances, estimated zone centroids were placed as 
close as possible to the centres of population gravity. In some 
cases the population gravity centroid varied markedly from the 
physical centre. The aim was to better reflect the actual distances 
travelled by people living in the zones to get to each of the 
centres. Travel distances were held fixed over the study period. 
In reality, changes would likely occur due to changes in road 
alignments, shifts in population centroids, and so forth. 


ie Income Estimates 


Separate estimates of per capita expenditures for Food and 
DSTM categories for each of the four study years were taken from 
the 1969 study. Income data appear to be the least refined of all 
data collected. Separate estimates had been made only for the 
four broad areal units constituting the total study area -- 
Belleville, Secondary Trade Area, Tertiary Trade Area A, and 
Tertiary Trade Area B. 


All 44 zones within Belleville were assigned a single average 
per capita figure for Food and for DSTM expenditures in line with 
the 1969 report, notwithstanding the variations in zonal per capita 
incomes that could be expected. 


Estimates for the remaining 13 zones also had to be aggregated 
for each of the three remaining areal units and each of the 13 zones 
assigned the corresponding average expenditure values. Projected 
per Capita Food and DSTM expenditures for 1976 and 1981 were 
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calculated in real dollar terms to avoid the complication of the 
effects of inflation on sales forecasts. 


Sa Food and DSTM Expenditures 

Estimating inflow or outflow of retail expenditures across the 
area boundary was difficult. It was estimated that 100 percent of 
the per capita expenditure on Food would likely remain within the 
study area. This assumption would not hold, however, for DSTM goods. 
The 1969 study estimated a 1971 outflow of DSTM spending from the 
four areal units of 20 percent for Belleville, 25 percent for the 
Secondary Trade Area and 30 percent for Tertiary Trade Areas A and 
B. Comparable levels for 1966 proved to be about five percent 
greater across the board. 


With construction of additional planned centres projected for 
1976 and 1981 it was assumed that the outflow could be reduced by 
a further five percent, but little hope was held for any further 
reduction. A certain amount of spending on specified classes of 
goods will probably always occur outside the Belleville study area 
by mail order and in nearby large metropolitan centres such as 
Toronto or even Kingston. 


9. Estimates of Centre Sales 


A further step was to determine what percentage of total 
potential expenditure was available to the centres being studied. 


Ratios of the total amount of floorspace for all the six or 
seven centres in the study to the total floorspace in the study 
area were calculated separately for Food and DSTM for each study 
year. 


The ratios were then employed to reduce per capita expenditures 
for Food and DSTM in each of the zones for 1966, 1971, 1976, and 
1981 to represent the actual amount available. 


When the total population in each zone was multiplied by its 
corresponding per capita expenditure, the 1966 estimate for Food 
and DSTM purchases was $49,500,000. The available floorspace 
multiplied by sales-per-square-foot ratios of $40 and $100 for DSTM 
and Food in all centres (except for one new planned centre where $60 
and $110 were applied) was $49,400,000. 
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For 1971, the comparable figures were $61,300,000 and 
$61,800,000 respectively. 


In both periods expenditure and sales appeared to be in 
reasonable balance, and the estimates employed were considered to 
be independently verified. 


CONCLUSIONS - DATA SOURCES 


Problems with data sources in Belleville were extensive. The 
availability of information from the 1969 study proved to be 
extremely helpful in some ways and created great difficulty in 
others. 


The problem of defining boundaries was much simplified as was 
the calculation of outflows and inflows of sales across area 
boundaries. 


The inventory of retail floorspace was excellent, but the lack 
of household information meant reliance on per capita expenditure 
data which was considered less suitable for purposes of calculating 
spending characteristics in a population. 


The difficulty of estimating actual sales for individual centres 
remains the “béte noire" of retail studies of any type. 


Some concern exists about the constraints imposed by having to 
allocate identical per capita expenditure figures over four large 
trade areas to 57 zones. Similarly, the number of zones concentrated 
in Belleville compared to the much smaller number covering the 
balance of the area is recognized as creating an added problem. 

While a large number of zones is preferable, the distribution of 
zones and the size of the population in each needs to be carefully 
balanced. 


The experience of the Board with the 1969 study and the 
Subsequent Ontario Municipal Board hearings (discussed in Shopping 
Centre Dectstons: Evaluatton Guides) appears to have made the Board 
conscious of the need for careful consideration of the retail 
component in planning studies. As a consequence, the Board was 
able to provide the study team with help in the interpretation of 
data for their area and in the preparation of testable options for 
future years that is not often available in other jurisdictions. 
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It can be concluded that the sources consulted provided only 
the minimum information needed to operate the model. Considerably 
greater refinement in detail would have been preferred. The 
existence of the trade area analysis done in 1969 reduced the 
amount of time required for data preparation. In addition, it 
provided invaluable independent verification of many of the estimates 
used. 


Options Tested 


Retail models are assumed to have the capability of predicting 
the impact of new proposed shopping centre sites on existing retail 
centres or the implications of a redistribution of floorspace among 
the centres -- for example, a gradual reduction of the proportion 
of space in the central business district and an increase in the 
amount of space in suburban shopping centres. Several options 
were tested for the projection year 1976 and 1981 in exploring this 
capability of the model. 


Decisions were made as to what options the study should test 
in advance of initiating any tests on the computer. The criteria 
used included considering which retail centres should be studied, 
what additional retail floorspace in new or expanded centres could 
be envisaged by 1976 or 1981, and what combination of centres and 
variation in floorspace seemed to make the most planning sense. 
The result was a proposal by the Board of the options to be tested 
(Table 3-2). 


Identical options were tested for the two versions of the 
retail model employed in the study. Table 3-2 presents a complete 
Summary of all the centres considered by the various options 
Suggested to the study team and the corresponding total floorspace. 
The location of each centre relative to the population zones within 
the study area is illustrated in Maps 3-3 and 3-4. 


The discussion of results is limited to an analysis of the 
years 1966, 1971, 1976 - Option 1 and 1981 - Option 6. Full results, 
however, were made available to Belleville for all options listed 
in Table 3-2, 


The base year 1966 was used to evaluate the sales levels of 
the Belleville C.B.D., Belleville Plaza, Towers, Trenton, Picton 
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and Napanee. By 1971, the number of centres had increased to seven 
following the opening of the Quinte Mall in 1970. The new centre 
provided an increase of 186,800 square feet. All the original 
centres remained static except for the Belleville C.B.D. A minor 
reduction occurred in floorspace from 344,800 to 342,000 square 
feet. 


The anticipated changes for 1976 - Option 1 were confined to 
increases in retail floorspace for the Belleville Plaza, 98,600 to 
168,600 square feet, and the Quinte Mall, 186,800 to 309,300 square 
feet. 


For 1981, seven options were decided upon and tested by the 
two retail gravity models. 1981 - Option 6 provides the best 
illustration of results. Nine retail centres were considered, 
added centres being the Trenton New Centre and the Leasehold 
Proposed Plans. A further minor reduction is shown for C.B.D. 
floorspace. The Board wished to establish the consequences of 
Significant growth and further dispersion of retail space by 1981. 


The options included in Table 3-2 represent only a few of the 
many possibilities that could have been tested by either of the two 
versions of the retail model. The selection of specific options by 
Belleville was constrained, among other reasons, by the amount of 
land available or desirable for commercial development, tentative 
retail proposals already before the Board, the ownership of certain 
lands, the owners' aspirations for commercial use, and the possibility 
Of some population growth in other communities within the retail 
sphere of influence of Belleville that might reduce the need for 
retail developments in Belleville itself. 


Calibration 


The Huff model was calibrated on 1966 data by varying the 
parameter value from 0.0 to 2.0. The best fit between actual and 
calculated sales appeared to occur at A = 0.88. At a value of 0.0, 
the model produced a near perfect fit, the “bogus calibration" fit. 
The sensitivity of the model to distance was eliminated and sales 
were distributed from zones to centres in strict proportion to 
size. A proper rank-order was maintained for calculated and actual 
Sales for all centres, further justifying the assumption of a 
lambda value of 0.88 (Table 3-16). 
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The Cullen model calibrated on the 1966 Belleville data by 
employing two parameters, alpha (a) and lambda (A). The best fit, 
where actual and calculated sales were most closely correlated, 
was found at a = 0.975 and dX = 0.036. 


The importance of maintaining a proper rank-order throughout 
the calibration process has already been discussed. The 
correspondence of rank-order can be seen by comparing the columns 
labelled SYN SALES and ACT SALES in Table 3-17. 


Predictive Capability 


After the calibration process was completed in such a way that 
the rank-ordering was correct and a high degree of fit was obtained 
for 1966, both versions were pushed forward in time by putting 1971 
data in place of 1966 data (Tables 3-18, 3-19). 


An analysis of the results indicates how well each version has 
done in predicting the 1971 retail situation as it was foreseen in 
1966, that is, their predictive Validity. 


For 1971 both the Huff and Cullen versions were able to predict 
the retail situation with considerable accuracy. In each case the 
rank-order of centres changed Slightly. For example, in the Huff 
model the rank-order was slightly altered between the first and 
sixth centres. The shift in rank-order was not considered enough 
to invalidate the predictive ability of the models. 


Impact of New Centres 


The results of the use of the Huff and Cullen versions to 
evaluate the consequences of the introduction of new centres at 
Specified locations is illustrated in Tables 3-18 and 3-19. 


It was suggested by the proponents of Quinte Mall during the 
Ontario Municipal Board hearings that preceded approval of the 
centre that the opening of the new mall would help or at least not 
harm Belleville C.B.D. 


If Table 3-18 is compared with Table 3-16 it will be seen that 
in the Huff version Belleville C.B.D., Centre 1, has been allocated 
sales in 1971 in excess of the 1966 level in spite of the 
introduction of Quinte Mall, Centre 7, which was allocated about 
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oro.e4 million in sales. 


In discussion with the Board in 1972, it was confirmed that 
the C.B.D. merchants did not appear to have suffered a loss in total 
dollar sales. This was not true for all merchants nor was it true 
for the relative shares of business being captured by Belleville 
feo.D. 


The Cullen version showed slightly different results with sales 
for the C.B.D. down about $1.0 million. At the same time it was 
reported that the trade area boundaries shown in the 1969 study had 
expanded because of the cumulative attraction of the new Quinte Mall 
and an attractive and well merchandised downtown. 


C.B.D. shopkeepers were said to be serving a Significant 
number of customers from beyond the boundaries of the Secondary 
Trade Area, particularly towards Peterborough. This would account 
Soretne ability of the C.B.D. to increase its sales. 


In addition, the extended boundaries would have increased the 
retail expenditure potential for the whole area and the amount of 
business that the model would have available to allocate to centres. 


The availability of consumer shopping survey data suggested 
for the Oshawa case study would have been useful in estimating the 
extent of the shift of trade area boundaries that appeared to occur 
in Belleville. A careful analysis of consumer Shopping trips is 
Once more underscored as an essential part of retail modelling. 


The usefulness of the two versions of the retail model in 
evaluating the impact of new centres appeared to be validated by 
the results in the Belleville area. Improvements in data and 
further experience with the operation of the models elsewhere 
Provide additional evidence that this is generally true. 


The Belleville experience suggests the advantages to be gained 
by maintaining a careful watch on results when new centres are 
Opened in major municipalities over a period of two to five years, 
On a year-by-year basis. 


It would appear that the Quinte Mall was built to the right 
Size and at the right time. It would also appear that the 
cumulative attraction of the new mall and the merchandising 
Strength of the C.B.D. prevented the serious loss of C.B.D. sales 
that has been experienced elsewhere. The Belleville Board 
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has confirmed the shortness of the recovery period in discussion 


with €:.B.D. merchants. 


The options tested for 1976 and 1981 in Belleville proved to 
yield excellent results. The amount of floorspace suggested 
matched the requirements almost exactly. It was possible, therefore, 
to do a detailed but tentative analysis of the results of the 


projections. 


A preference was given to 1976 - Option 1 over 1976 - Option 2 
because of the better fit obtained between calculated sales and 
actual sales. The magnitude and distribution of retail floorspace 
used in 1976 - Option 1 better satisfied the projected demand 
resulting from income and population increases. A review of Table 
3-2 will identify the floorspace variations considered in 1976. 


For the Huff version, 1976 - Option 1 was tested using the 
calibrated lambda value of 0.88 and an actual sales-per-square- 
foot figure of $63. In two instances the desired rank-order 
between calculated sales and actual sales shifted between 1971 and 
1976 (Table 3-20). 


The projected expansion of the Quinte Mall and the Belleville 
Plaza would appear to satisfy the majority of the potential retail 
demand. In addition, a satisfactory sales level is maintained 
among all seven centres. There would seem to be no justification 
for any additional large-scale shopping facilities before 1976, 
assuming that the anticipated expansion occurs at the two selected 
sites and that population and incomes increase as predicted. 


A review of customer distribution and calculated sales among 
the seven centres reveals an interesting trend. As suggested in 
1971, Belleville C.B.D. will capture an even smaller percentage of 
the total consumers (30.69%) and of calculated sales (31.72%). 
Meanwhile, the floorspace additions to the Belleville Plaza and 
Quinte Mall would increase their percentage share of consumers 
from 10.16 to 14.30, and 16.90 to 23.80 respectively. Similar 
increases would also occur in the calculated sales distribution. 


In 1976 - Option 1, further expansion of the Belleville Plaza 
and Quinte Mall could be justified on the basis of the projected 
income and population forecasts used in this study with some degree 
of certainty of a non-deleterious impact on the existing retail 
structure. Any further large-scale shopping centre proposals would 
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have to be considered premature in light of the 1976 - Option 1 
results. 


The summary for the Cullen version for 19/6. 0p tion 14s 
shown in Table 3-21. By the end of the five-year period total sales 
increased by approximately $14.5 million, caused by an increase in 
income and population. The fit for this option is again quite good 
but not as good as that obtained in the previous two periods. 


The increase in demand for retail floorspace was met by two 
changes in existing retail facilities. The Belleville Plaza 
increased its floorspace by 70,000 Square feet making it the third 
largest centre. The size of the Quinte Mall increased by 122,500 
Square feet to a total of 309,300. Quinte Mall at this period was 
approaching the size of the C.B.D. 


It is difficult to determine the effect of a Single centre 
increase when more than one centre is Changed within the same time 
period. If individual effects are to be determined, they must be 
examined in separate options. In this Study the effects are 
considered in terms of both increases occurring simultaneously. 


The C.B.D.'s decline in importance is reflected in the customer 
and cash flow matrices (Tables 3-26, 3-27). C.B.D. ranked first 
meves9 of the 57 zones -- a decrease of 12 from 19712, “ine theszones 
where the C.B.D. was displaced in rank, its position was relinquished 
almost exclusively to the Quinte Mall. 


It is apparent that when adequate information is available the 
retail model shows considerable potential for examining the probable 
consequences of the introduction of new centres at stated years in 
the future. The analysis of Belleville's 1976 - Option 1 appears 
to provide a convincing illustration. 


Projections to 1981 were also attempted. Seven options were 
actually tested. The best fit was produced by 1981 - Option 6 
which is discussed here. 


Using the Huff version it would appear that by 1981 there is 
Substantial indication that an Opportunity will exist for the 
further development of large-scale shopping facilities (Table 3-22). 


The summary site information reveals that two plazas -- Trenton 
New Centre of 100,000 Square feet and Leasehold Proposed Plaza of 
121,400 square feet -- could be justified and operating by 198] 
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without any detrimental impact on the sales levels presently 
enjoyed by the seven existing centres. 


The level of sales for the existing centres would increase, 
maintaining an acceptable sales-per-square-foot ratio. The 
distribution of consumers and total calculated sales would decrease, 
however, aS a consequence of increasing the number of centres. In 
effect, the percentage share of the centres would decrease because 
customers and sales would be redistributed among nine instead of 


seven centres. 


According to the predictions of the Huff version, 1981 - Option 
6 seems to maximize overall projected customer shopping satisfaction 
and best provides a combination of centres to satisfy both regional 
shopping patterns and intra-urban shopping trips. Any other changes 
in the retail structure in the interval before 1981 which have not 
been included in the options tested would obviously invalidate this 
conclusion. 


The testing of 1981 options for the Cullen version confirmed 
that 1981 - Option 6 appeared to best satisfy the increased retail 
Sales demand of $12 million which was generated in the previous 
five years by population and income increases (Table 3-23). The fit 
is again quite good but not as good as in the three previous 
analyses. 


Two centres have been added to the existing seven in order to 
meet projected floorspace needs. Trenton New Centre has a floorspace 
of 100,000 square feet which places it sixth in terms of size. The 
Leasehold Proposed Plaza ranks fifth in size with a floorspace of 
121,400 square feet. 


It is again necessary to examine the impact of two simultaneous 
changes in retail facilities. The two new centres would appear to 
be able to obtain 15.4 percent of the retail market with seven 
percent going to the Trenton New Centre and 8.4 percent to the 
larger Leasehold Proposed Plaza. 


All of the original seven centres were shown as suffering a 
decline in percentage share of the retail market. The most 
significant losses were allocated to the Belleville C.B.D., a 3.7 
percent drop, Trenton Centre, a four percent decrease, and the 
Quinte Mall which declined 3.6 percent. The Napanee Centre seems 
least affected with a decrease of only 0.5 percent. 
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In terms of sales per square foot, both of the new centres 
are projected to come in with below the overall] average of $63.50 
per square foot. The Leasehold Proposed Plaza reaches only $60.60 
per square foot. 


The most significant change is a projected decrease of $12.6] 
per square foot suffered by the Trenton Centre which seems to be 
caused by the Trenton New Centre. This Suggests that a possible 
Over-expansion of retail floorspace in Trenton could result by 1981. 


The two largest centres, the Belleville C.B.D. and the Quinte 
Mall, remain quite stable, each showing a $0.10 increase. Napanee 
Showed an increase of $5.88 which is Significant in that it now 
has the greatest sales-per-square-foot figure -- $80 per square 
foot. At this point Napanee may be understored in relation to the 
Overstored projection made for Trenton. Thus, there are indications 
that Napanee could stand on expansion of its retail space by 1981. 


The ranking of the Belleville C.B.D. in customer and cash flow 
matrices remained constant in its number one PeSTEDONM i N39: Od tthe 
97 zones. The Quinte Mall was again most commonly found in the 
number two position and, when the C.B.D. slipped from its number 
One ranking, it was usually replaced by the Quinte Mall. 


In the years projected by the Cullen version of the retail 
model, it appears that the dominance of Belleville's Cb. can 
Still be maintained. The proposed new centres and the expansion 
of existing centres foreseen by the Board would not leave 
any centre in an untenable competitive situation although Napanee 
appears to be short of space and Trenton oversupplied. 


The outputs of the model contain income and population increases 
which determine consumer demand. To maintain the validity of any 
Prognostications of this type, continued monitoring of current 
Changes would be essential. 


Levels of Service 


The Belleville study provides a clear indication of the 
usefulness of the measures of level of service provided by the 
model, 


In the discussion of 1981 - Option 6 it is apparent that the 
allocation of retail floorspace to Trenton and Napanee would result 
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in an overstored and an understored condition respectively. 


In this instance the customers in Napanee would not be 
receiving a level of service comparable to what they had previously 
experienced. The use of the retail model facilitates this type of 


determination. 


The use of the average-trip-length measure of level of service 
was employed extensively in the Belleville analysis. The options 
for 1976 and 1981 all produced enough total retail floorspace to 
match anticipated demand. 


It was possible to compare the average trip length produced by 
the distribution of retail floorspace reflected in each option to 
the average trip length produced by all other options. This measure 
was employed by the study team as one basis for selecting an option 
to use for purposes of illustration -- the one with the best overall 


performance. 


Comparative Measures of Retail Productivity 


The Belleville study produced numerous illustrations of the 
usefulness of the measure of retail productivity in analyzing 
particular options or comparing one option against another. This 
ability of the model would appear helpful to the retail analyst 
working with a shopping centre developer or planning board. 


Customer and Expenditure Flows 


The use of zonal expenditure information and zone-to-centre 
flows of customers and expenditures has already been illustrated in 
earlier discussion of the Belleville applications. A sample of a 
Huff neighbourhood summary is included as Table 3-24. The tables 
produced by the Cullen version of the model are shown as Tables 3-25, 
3-26, 3-27. The matrices provided with the Cullen results are easy 
to follow and have proved to be most helpful when trying to trace 
faults in the summary tables or to ensure that improbable shopping 
trips have been adequately constrained in the calibrated version. 


Conclusions - Belleville 


The sources of data in the Belleville study have already been 
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discussed. The Belleville case study provided an Opportunity to 
examine the information produced by the model in considerably more 
detail than was possible in Oshawa. 


The calibration procedure worked out easily for both versions 
and the test of predictive Capability showed enough precision that 
short-term (five year) forecasts could be made with some confidence. 


The options suggested for testing in Belleville proved to match 
the future retail floorspace requirements in all instances. The 
four options chosen for discussion from among the eleven actually 
tested show quite clearly that the retail model is capable of 
giving a perceptive indication of the past and present retail 
situation in Belleville. 


Moreover, it was illustrated how planning decisions about the 
Provision of retail facilities for the future can employ estimates 
of consumer travel as a means of ensuring adequate levels of 
customer service and centre Viaba lity. 


Measures of retail productivity also proved to be useful 
indications of customer service within the limits of currently 
available forecasting methods. 


Future changes in population, consumer travel patterns, per 
capita income, per capita retail expenditures, attractiveness 
mactors etc., all have a profound effect on the suitability of 
retail locations and shopping satisfaction. 


The application of the retail model in Belleville appears to 
Show quite adequately how it can assist in reducing uncertainty 
in retail planning decisions and to suggest quite strongly that the 
retail model does have a role to play in ensuring a responsive and 
vital retail segment in the community. 


CONCLUSIONS - CASE STUDY MUNICIPALITIES 


Sources of Data 


1. Floorspace Estimates 
In both Oshawa and Belleville it was confirmed that adequate 


data sources were available to operate both versions of the retail 
model tested. An up-to-date retail floorspace inventory was 
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available in both instances, and this proved to be the single most 
important set of required data. Several tests of the reliability 
of the inventories were carried out including field survey checks 
of a few centres and confirmation that sales per square foot for 
selected trade categories were within acceptable bounds. 


he Zones & Boundaries 


The selection of zonal boundaries appears to require some care. 
In both instances the boundaries of the study area were largely 
predetermined by previous retail studies. In the case of Oshawa 
the zones were fairly evenly distributed throughout the area. In 
Belleville the concentration of zones within the city created an 
even population distribution throughout the trade area, but a less 
than satisfactory spatial distribution. 


ae Population, Income and Expenditure Estimates 


Zonal household figures and expenditures in Oshawa were estimated 
with a high degree of reliability. Household information was not 
readily available in Belleville on a zonal basis. A compromise 
solution was adopted involving the use of zonal population figures. 
Per capita expenditure figures could not be disaggregated zonally 
in the Belleville study. The Oshawa data would appear to be 
superior to the Belleville data in these respects. 


4. Distance Estimates 


Distance estimates for Oshawa were available from Ministry of 
Transportation and Communications survey data. For Belleville, 
distance estimates were scaled directly from maps, a less reliable 
procedure. A test of the straight-line distance option in the Huff 
model proved difficult to evaluate in relation to travel-time matrix 
measures. However, there seemed to be little apparent difference in 
the allocations of sales in Oshawa when employing measures of 
straight-line distance calculated internally by the Huff program 
from grid co-ordinates. 


el Centre Sales 


There were no significant differences in the methods used to 
approximate actual centre sales between the two case study 
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municipalities. The averaging method used in the Huff version 
tested would appear to be less satisfactory than the variation in 
centre sales per square foot available in the Cullen version. This 
discrepancy has since been corrected. The approximation of actual 
centre sales still remains an area in need of further Study... his 
is particularly true when the potential sales of contemplated new 
centres are being evaluated. 


6. Options Tested and the Impact of Future Floorspace 


A comparison of the two municipalities reveals the significance 
of the ability to estimate future floorspace requirements in relation 
to future distributions of population and retail centres. The 
Belleville table of options allowed a detailed comparison and 
evaluation of alternatives. In the case of Oshawa, the potential 
of the model as a means of projecting future floorspace requirements 
was demonstrated. 


Tie Comparative Measures and Flows 


In Belleville the opportunity to test the comparative measures 
of level of service, productivity and customer and expenditure 
flows suggested how the measures could be usefully employed in 
evaluating alternative locations and centre sizes. In both 
municipalities it was apparent that the customer and expenditure 
flow matrices of the calibrated version could have been checked 
against surveys of customers' shopping behaviour, had they been 
available. The Huff version of the model allows for the 
calibration of the model against consumer surveys to ensure that 
the resulting solution matches observed effects. 


8. Comments 


The comparative testing that was conducted appears to confirm 
the principal claims of the usefulness of the retail models. The 
results reproduced in the case studies have been discussed in some 
detail. 


It is obvious that it would be worthwhile to repeat the studies 
to make a number of significant adjustments in the data that were 
Provided to the model, and to introduce further variations by a 
full testing of each of the options available within each version. 
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The usefulness of doing so was not contemplated within the adopted 
research design. 

In any case, subsequent operational applications of the model 
in places such as Oxford County, Regional Waterloo, Regional 
Niagara, Kingston and London benefitted greatly from both the 
advantages and shortcomings that were revealed in the applications 
in the case study municipalities. 
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GLOSSARY OF TERMS - TABLES 3-3 THROUGH 3-27 


SITE The number of the centre in the sequence 
in which they appear in the data provided 
to the model. 


SIZE The size of the centre in Square feet of 
gross leasable area or its equivalent for 
Food and DSTM merchandise. 


FOOD The square footage of Food floorspace in 
the centre. 


DSTM The square footage of DSTM floorspace in 
the centre. 


AVTL The average trip length for customers 
allocated to centres reflecting trip 
distances weighted by their associated 
customer flows. The average trip length 
over the entire study area is shown as a 
total. 


XS,YS The X and Y grid coordinates of each 
centre measured in decimal inches from a 
predetermined origin. 


CUST The number of customers allocated to each 
centre on the basis of the probabilities 
calculated by the program. 


eat SLS, The calculated or synthesized sales. The 
SYN SALES 
product of the number of customers from 
each zone and their average expenditures 
On Food and DSTM goods. 


ACT SLS The approximation of actual sales derived 
as a product of floorspace and sales per 


Square foot. 


CS/ Fi 


AS/FT 


HH SH 


CA$SH 


AC$SH 


GLOSSARY OF TERMS (CONTINUED) 


The calculated sales per square foot, 
derived by dividing SIZE by CAL SLS. 


The actual sales per square foot derived 
DY diViding SIZE DYeACI eso. 


The share of households or populations 
attracted to the centre from the total 
households or population in the area. 


The estimate of the share of calculated 
sales attracted to each centre. 


The estimate of the share of actual sales 
attracted to each centre. 
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4 General Conclusions and 
Areas for Further Research 


The applications in the two case-study municipalities have 
demonstrated the power of the two versions of the retail model to 
help the urban and regional planner arrive at public policies on 
shopping centres and at shopping centre decisions. 


The Cullen version of the model appears to have somewhat 
‘Superior operating characteristics. The predictive capability of 
both versions of the retail model appears to be at least equal to 
or better than that obtained by conventional methods. 


Because retail models and trade area analysis rely on similar 
data bases and face the same problems of estimation, the demonstrated 
predictive ability of the retail model gives it a powerful advantage. 


In effect, the research has confirmed the postulations that 
retail models can evaluate a wide range of alternative arrangements 
for providing retail shopping facilities under varying assumptions 
With a high degree of reliability. Generally, data is available 
to operate the models and perceptive measures of forecasted results 
are obtained. 


Other less dramatic but important conclusions about retail 
models were reached. 


a. The problems of establishing a sound data base are quite clear. 
Any municipal or regional planning study is only as sound as 
the quality of the information on which it is based, and the 
application of the retail model is no exception. 


High-quality data for past and current years are required 
whether a site evaluation of the step-by-step type is to be 


ao eee 


employed or a retail model is to be used. 


Subsequent applications of the retail model in real-world 
planning studies in Ontario have substantiated the notion that 
it is essential to do a trade area analysis as a preliminary 
step in approaching a model application. A survey of consumer 
shopping is particularly important in this regard. 


Household Incomes. Finding sources of income information for 
small statistical areas was a serious problem which was only 
partially overcome in both case study areas, although greater 
success was achieved in Oshawa than in Belleville. 


Income information will improve only in part once 1971 census 
information from Statistics Canada becomes available. The 1971 
census will provide average household income information for 
statistical areas for the principal urban centres based on a 
One-third sample. Estimates for rural areas will remain 
difficult. Moreover, it will still be necessary to derive 
estimates for intercensal years. 


Other sources may have to be employed. These include the 
personal knowledge of local and provincial planning staffs, 
taxation statistics, and the co-operation of other bodies that 
possess some income information such as local credit bureaus. 


Periodic random sampling could provide continuous monitoring 
of income for small statistical areas if the benefits of the 
analysis could be shown to offset the costs. If the income 
data were collected as part of normal Survey work the added 
cost would be very small. 


Centre Sales. With rare exceptions, annual sales for individual 
shopping centres are presently not available from any published 
source. The disclosure rules of Statistics Canada quite 
properly preclude their release. Users of the retail model 

are in no worse position than trade area analysts in this 
respect. A ratio of floorspace to sales was used in the case 
studies as the best substitute for accurate information on the 
annual sales of existing centres. 


Elsewhere, planning staffs have been able to obtain independent 
confirmation of their estimates through local retailers, annual 
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reports of shopping centre operators and other financial 
institutions. 


The analysis of the multiplicity of options tested in the two 
case-study municipalities amply demonstrates the ability of 
the model to handle a very large number of alternative 
patternings of retail distributions simultaneously in a short 
time at low cost. 


The amount of computer time required to produce a summary site 
table of the type used for illustrative purposes in the case- 
study chapters was only a matter of a few seconds. For example, 
the computer time required to run the Huff version was less 

than five seconds. 


All the tables produced in analyzing the options tested for 
Belleville or Oshawa using the Huff program amounted to only 
minutes. The full version of one test for one option for the 
Cullen program took only 109 seconds at an estimated cost for 
computer time of $1.42. 


The principal amount of time and effort required to establish 
an adequate high quality data base is measured in terms of 
weeks or even man-months. The same statement can be made 
about alternative patternings of residential areas, the road 
network, households, income and so forth. 


The model demonstrated its usefulness in determining whether 
adequate provision has been made in a plan for future retail 
requirements. The Belleville options achieved equilibrium 
solutions in each of the study years. In Oshawa, the results 
of the options tested provided a clear picture to the planners 
of the direction to be followed in establishing an adequate 
supply of potential sites. 


The impact of the introduction of added retail space at a 
given location on the remainder of the shopping centres is 
clearly shown in terms of customer flows, cash flows, sales- 
per-square-foot ratios and market shares. 


Some measure of the degree of acceptability of one site over a 
number of alternative sites -- where more than one site is 
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under consideration and there appears to be justification for 
only a single centre -- can be obtained by examining the 
relative size of the average trip lengths. It can be assumed 
that the pattern of shopping centres with the lowest average 
trip length contributes most to customer satisfaction by 
providing convenient access to shopping facilities. 


oe Where two proposed centres are located side-by-side, and there 
is ample justification for only one centre, criteria beyond the 
value judgements introduced by trade area analysis or the retail 
model must be introduced by local or regional government decision 
makers in any case. 


fe The use of either version of the retail model obviously assists 
in establishing the initial desirable size and future phasing 
of expansion for a single centre or for a complex of centres. 
Problems of overstoring or understoring in one locality are 
quickly revealed. 


8. The research program undertaken by the study team generally 
confirms the validity of suggestions for further research that 
have been put forward by others. As well, priorities for 
further specific areas of research have been considerably 
clarified. 


AREAS FOR FURTHER STUDY 


First priority in further research should be given to 
improvements of a practical, and more immediate, nature aimed at 
further modifying the two Programs used in this study for planning 
purposes. 


Second priority could then be given to research of a theoretical 
and behavioural nature to strengthen the underlying conceptual basis 
of the retail model and to incorporate more recent thinking of 
other researchers. 


Operational Improvements 


A number of substantial improvements and refinements have been 
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incorporated into the two programs presented in this report. In 
each case, the program output has been considerably expanded, and 
this is particularly true of the Cullen version. 


Moreover, the outputs for each program are now more readily 
comparable, so that a much clearer and more detailed picture of the 
actual situation for any retail centre, existing or projected, is 
presented. Interpretation of the results for policy purposes is 
greatly facilitated. 


There are a number of interesting possibilities for further 
modification. The following refinements, for example, could prove 
to be valuable. They could not be attempted in this study due to 
time constraints. 


ie The incorporation of a neighbourhood printout in the Cullen 
version. This is one important advantage the revised version 
of the Huff and Blue program possesses. The customer and 
cash flow matrices given in the Cullen program are of great 
value. Nevertheless, the form of data output in the neighbour- 
hood printout available in the Huff and Blue version is easier 
to read for some purposes. 


(ag The modification of the two programs to use identical data 
decks if further comparative testing of the usefulness of each 
is to be continued. At present, each program requires virtually 
the same data but arranged somewhat differently. Such an 
adaption would result in a considerable saving in time and 
keypunching costs, especially when a large number of zones and 
centres are involved. Some effort has already been devoted to 
this area. The remaining changes would likely require some 
Program restructuring. 


3). A review of calibration procedures and the possibility of 
introducing another calibration process into the Huff and the 
Cullen versions. 


4. The incorporation of further refinements to include all levels 
in the hierarchy of centres and all expenditures suggested in 
recent advances made in the U.K. in multi-centred hierarchial 
versions. 
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Introducing operational modifications would help the user 
without altering the underlying concepts in any way. 


On the other hand, the existing programs already allow a 
reasonably sophisticated analysis to be made. 


The more reliable the data base used for analysis, the more 
reliable the two programs become as planning tools. 


Thus the argument for making further program revisions is 
perhaps not as strong at this time as the case for incorporating new 
theoretical foundations into the retail gravity concept itself. 


Theoretical Improvements 


The retail trade area analysis approach has many underlying 
weaknesses as a working model. It is oriented almost entirely to 
the single centre study as the basic unit of analysis. It is of 
little value for describing either long-range or multi-centre retail 
strategies. Analytically, it is entirely dependent on the use of 
directly observable phenomena such as accessibility and traffic flow, 
population and its distribution, income, and existing and potential 
competition for describing the viability of any site or the probable 
spatial configurations of retail spending. No acknowledgement is 
made of less directly measurable, more behavioural elements which 
might affect a consumer's choice in a significant way. 


The methodology is satisfactory from the viewpoint of a 
developer interested in a reliable estimate about the soundness of a 
particular retail development proposal. It serves as a guideline for 
deciding if a selected location is appropriate as an investment. 


As a planning tool it is much less satisfactory. The step-by- 
step procedure requires the use of privileged information and a 
series of shrewd estimates at each stage. It is more characteristic 
of an art than a soundly based scientific approach. 


Retail models are also static, and therefore do not formally 
predict future shopping patterns or capacities. They do provide, 
however, a quantitative analytic base from which predictions may 
be obtained by projecting certain variables. 


It is possible to adjust individual parameters and variables 
within the models to discover the possible effects of time and to 
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test and generate alternative plans. This flexibility is a very 
attractive characteristic of the more recent models and makes them 
much more valuable as planning tools than trade area analysis. 


The gravity model concept is not without weaknesses. It is 
heavily dependent on the same variables as trade area analysis. It 
has been suggested that success in formulating a model of individual 
Or mass behaviour requires the inclusion of variables other than 
strictly economic factors. However, to do so would necessitate a 
much sounder understanding of individual choice behaviour than now 
exists. 


The consumer does not seem to follow a wholly rational process 
when travel costs are traded off against the shopping centre's 
attractiveness. Habit and inertia have been suggested as being 
among the dominant influences. 


Current models of analysing the shopping process have still 
only begun to approach the question of individual consumer behaviour. 
Macro-scale models for forecasting or planning purposes have great 
potential but to date are still only rough tools in the absence of 
knowledge of the true underlying variables. 


Although encouraging results have been obtained, current 
versions may be little more than historical accidents in view of 
the absence of a theoretical connection between their empirically 
determined weights and exponents and the corresponding behavioural 
variables on which they rest. 


Intervening Opportunities Model 


Hayes discusses a further form of retail model proposed by 
Harris (Hayes 1968). This too is a probabilistic model. Its main 
postulates are that the number of shopping trips or sales from an 
Origin in zone i to a destination at shopping centre j are 
proportional to the number of intervening opportunities. 


Unlike the gravity model, interzonal distance does not appear 
explicitly. Instead, possible destinations are ranked in order of 
increasing impedance from the origin. 


Although the model in practice resembles the retail model that 
has been tested here, some interesting improvements appear possible 
Once further refinements and testing are done. 
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lg The model contains a specific provision to vary demand for 
different population classes. 


(ap The model contains a provision to vary consumer behaviour in 
accordance with the density of shopping opportunities. 


Ss Two kinds of market are recognized, home-based and work-based 
Origins. 
4, Wider description of consumer preference is possible. 


ae Seven parameters are contained within the model although data 
for calibration purposes are not usually available and, 
therefore, when simplified to an operational form, the reduced 
version approximates a gravity type model. 


Future Directions in Retail Modelling 


Retail models, individually or cumulatively, provide much 
information to serve as guidelines towards decisions on retailing. 
There is still ample room, of course, both theoretically and 
practically, to improve on present methodologies. 


The results of preliminary attempts at microanalytic modelling 
of consumer preferences and shopping patterns are only now beginning 
to appear. The work of MacKay is instructive in this regard (MacKay 
1972.6 oy 30s 


An underlying weakness is the lack of behavioural information 
to support the variables and parameters used in gravity models. Two 
particularly unsatisfactory gaps in the information available for 
retail model building include: 


12 a need for a better understanding as to how shopping will be 
carried out in the future, and 


ce a need for further knowledge about the shopping behaviour of 
individuals and small groups instead of relying on the actions 
of large aggregates of people. 


le Changes in Retailing 

A forecast of future shopping activity would have to try and 
answer two important questions. First, how will changes in the 
retail and service trades affect the "attractiveness" of different 
Sizes and types of shopping centres? And second, how will changes 
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in the way that shopping is carried out be reflected in the 
consumer's trade-off between the cost of a journey and the purpose 
of the journey. 


2. Consumer Behaviour 


A better understanding of how people really perceive their 
choices in shopping has been advocated (Thompson 1966). 


Thompson has concluded that a more interdisciplinary approach 
is needed in attempting to identify promising new avenues of 
theoretical and empirical investigation using modelling techniques. 


It must be recognized that the geographic distribution of 
retail sales is the result of the summed reactions of many 
individuals, each free to make decisions given the perceived and 
real constraints with which the consumer is faced. 


To build such notions into retail models would in effect call 
for a halt to further immediate research along the traditionally, 
empirically-oriented areas of investigation which the present 
retail models have followed. 


A closer link would be needed between the empiricists and 
theorists to establish a truly inter-disciplinary viewpoint. 
Thompson suggests that this can be achieved if the following 
hypothesis is adopted: 


The fundamental factor affecting the geographic 
distribution of retail trade is the manner in 

which consumers organize their perceptions of 

the internal environment with which they are faced. 


This would call for less emphasis on census and other published 
data sources as the principal basis for the operation of explanatory 
models, and more emphasis on consumer shopping strategies and 
Survey research. 


Studying how individuals interpret and make decisions about 
complex phenomena in experimental situations may well be of 
potential value in formulating more meaningful theories of shopping 
behaviour. This in turn, may serve as the basis of further attempts 
at model building and quantification in the area of retail area 
research. Thompson cites several examples to illustrate the 
importance of this approach including the work of Lynch on the 
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perception of the urban environment. 


If the orientation of the individual is demonstrably 
subjective, this could be important to the problem of analysing 
shopping behaviour and the geographic distribution of retail 
patronage. Much could be gained from attempts to explain perception 
differences in terms of such measurable consumer attributes as 
length of residence, education, income, age, sex and social class. 


Summary 


It may be that advances in retail modelling can only be 
achieved by integrating knowledge from many fields and incorporating 
this into the model using the statistical mechanical methods of 
entropy. 


This calls for the development of a proven theoretical basis 
to support the variables and parameters used in the application of 
quantitative analysis -- both on a micro and macro level -- when 
applying modelling techniques. 


Such techniques would likely be used with even greater 
confidence if a multi-disciplinary approach could be devised for 
future studies of retail location and consumer behaviour. 
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Steps in Using the Programs 


The purpose of the technical appendix is to outline the two 
computer programs used in the study as a means of familiarizing the 
reader with the operating procedures involved. 


Obviously, the technical aspects of computer applications can 
appear somewhat confusing to the reader who is not acquainted with 
programming languages and computer techniques. 


While the appendix has been included largely to assist those 
interested in actually conducting similar studies, the format 
adopted should not prevent those who must interpret the results 
from examining the procedures. 


A full computer program consists of three parts. 


First, there is the information that the computer must be given 
in advance so that it will have something with which to operate. 
This is called the input data. 


Second, there is the set of instructions to the computer on 
What is to be done with the input data when it gets them. This is 
called the source program. 


Finally, there is the printed result of the computer operations. 
This is called the output. 


Input data consist of all the variables that are to be considered, 
_the parameters of the model, the options that are available to the 
user and similar kinds of information. 


In the case of the retail model, the input variables comprise 
' Such things as the number of centres under study, centre sizes, 
distances from residential areas to retail centres and retail 
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expenditure information. Output variables, on the other hand, are 
the results of the computations made by the computer. 


Both input and output variables must be carefully specified to 
ensure that the results are reliable. If the information required 
for the program is not properly interpreted, the results could be 
useless. Similarly, results that are otherwise reliable could be 
seriously misconstrued if output variables are misread. 


THE HUFF AND BLUE PROGRAM 


The program developed by David L. Huff and Larry Blue has been 
revised in the course of the study and used to test the Huff version 
of the retail model. 


The copy of the source program is accompanied by a simple set 
of step-by-step explanations serving as a basic guide to the 
prospective user. 


The variables in the model are described and the names assigned 
to them for programming purposes are identified. 


Here are detailed descriptions of the input data requirements 
and all the variables in the program's output. 


The Source Program 


A copy of the source program will be found in Annexure A. An 
effort has been made to reduce the program listing to a form 
which can be more readily interpreted by the reader. 


Each step is accompanied by explanatory notes about the nature 
of the operation being performed. A simplified flow diagram has 
been added in the margin to illustrate the sequence of operations. 


The source program is written in FORTRAN IV language and is 
Suitable for use with any WATFOR compiler. Modifications may be 
necessary if it is run on another system. 


To make the program somewhat easier to follow, si tis® use fu od 
compare the basic Huff formulation with its corresponding translation 
into program language. For instance the notational version of the 
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model 


appears in the program version as 


The mathematical notations have been changed to machine readable 
variable names to match the conventions of the programming language 
being used. Explanations of each of the variables employed in 
translating from notational conventions to program conventions are 
found in Figure A-1. This will serve as an introduction to the way 
in which variables throughout the program are named and defined. 


Description of the Input Data 


The following punch cards are required as inputs into the 
Program. The cards are described in the sequence in which they must 
appear. A number of options are available to the user and their use is 
described in some detail at the point at which the option can be invoked. 


Title Card 
Column Notation Definition 
1-72 TUPLE As many as 72 alphanumeric characters may be 


employed to identify the study being undertaken. 
Any combination of characters, letters or 
numerals, up to a maximum of 72, may be used. 
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Control Card # 1] 


Column Notation 


1-10 ENI 
11-20 ENJ 
21-30 EN] 


31-40 BRS 


41-50 ER SE 


Definition 


A numerical value corresponding to the number of 
zones or statistical areas within the study 
region. 


A numerical value corresponding to the total 
number of retail locations within the area 
under study in any given run of the program. 
This value cannot exceed the value assigned 
to J in the DIMENSION statement. 


A numerical value corresponding to the maximum 
number of iterations specified for the lambda 
search procedure. A suggested safe number to 
use is 50. On the other hand, if a value for 
lambda is specified and no parameter search 
procedure is desired EN] must be set equal to 
0. Consumer shopping survey data are required 
if the parameter search procedure is to be 
initiated. 


A value representing the minimum lambda 
increment employed in the Fibonacci search 
procedure. The value suggested by Huff and 
Blue is .00001. If no search is being 
conducted a 0 should be specified. 


The minimum error increment to be used in the 
Fibonacci search procedure. Huff and Blue 
Suggest a value of .0000001. Again, a 0 
Should be specified if no search is to be 
executed. It should be noted that not all 
computer systems may be able to accept values 
for EPSL and EPSE as small as those suggested 
by Huff and Blue. Where this is the case and 
a search procedure is desired, a reduction of 
several decimal places should overcome this 
ditiiaLculeys 
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NOTATIONAL AND PROGRAM VARIABLES - HUFF 


Mathematical Program 
Notation Notation Definition 


The probability of a consumer or 
household from a zone i shopping at 
a retail centre j. 


The number of consumers or households 
ined ZONeG i. 


The average annual expenditure per 
consumer or household budgeted for 
a given class of goods in a zone i. 


The total retail expenditure 
available in a zone i. 


The sales attracted to a centre j 
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from a zone i. 


A measure of the attractiveness of a 
centre j, in this case floorspace. 


A measure of the travel distance, time 
or cost from a zone i to a centre j. 


The number of retail centres in a 
given study area. 


A parameter estimated empirically to 
reflect the effect of travel distance, 
time or cost on various kinds of 
shopping trips. 


Fig. A-1: Notational and Program Variables - Huff 
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RATE 


RFIND 


Control Card # 2 


Column 


1-10 


11-20 


Notation 


SITEN 


AVIND 


A numerical value corresponding to some 
predetermined straight-line distance 


conversion factor, such 
etc. For example, RATE 
average dollar cost per 
applied to all shopping 
conversion is desired a 
is simplest to set RATE 
distance data are to be 


as time, miles, cost, 
could represent an 
unit distance to be 

If no such 

1 must be punched. It 
equal to | if actual 


trips. 


used, though the value 


of this indicator is not critical in this case. 


An indicator that, when assigned a value of 0, 
will bring about a calculation and printout of 
A O must 
be entered if the lambda search is to be 


R, the coefficient of correlation. 


executed. 
cause this step to be omitted. 
to enter a value of 1 only when no parameter 


Specification of the value 1 will 
It is necessary 


search is desired. 


Definition 


An indicator that identifies any retail site j 
to be increased in size incrementally. This 
indicator should be punched as an integer 

value corresponding to the position of such a 
centre in the array of location cards. For 
example, if the third centre listed was to be 
incremented then 3 would be punched in this 

A 0 must be specified if no centres are 
to be incremented. If SITEN is not equal to 0, 
ASC, a control indicator on the following card, 
may be left blank. However, if SITEN is 0, 

ASC must be specified. 


space. 


An indicator that, when assigned the value of 
1, brings about a recalculated neighbourhood 
and summary site printout using a determined 
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21-30 


31-40 


41-50 


PIT 


PNT 


CTRI 


average lambda over all zones. In this case 
the optimum parameter values for each zone 
derived by the search procedure are summed and 
divided by the number of zones to arrive at 

an overall average. If AVIND is assigned a 
value of 1, PNT, an indicator on the same 
card, should also be 1. AVIND must equal 0 if 
this average lambda option is not desired. 


An indicator that, when assigned the value uA 
brings about a printout by zone of the 
Fibonacci search for the parameter lambda. No 
such printout is given if a 0 is punched. If 
AVIND is equal to 1, PIT should be assigned a 
Valuerof 02 


An indicator that, when assigned the value Le 
results in a printout of the sales and customer 
distribution to each centre by zone. Such a 
printout can be very useful but is likely to 

be very lengthy and thus most costly. AO must 
be specified if such output is to be suppressed 
and only summary site information is desired. 
However, if the parameter search option is 
incorporated PNT should equal 1 as the zonal 
printout cannot be properly suppressed in this 
case. Similarly, if AVIND is assigned a value 
of 1, PNT should also be 1. 


An indicator that determines how many sets of 
Zone cards are to be read. There can be a 
maximum of 3 such sets. The first set must 
always be read and represents summary 
information for each zone. The second set 
represents actual consumer survey data and 
indicates the number of customers travelling 
from each zone to each centre. The third set 
represents actual distance data in terms of 
cost, travel time, miles, etc., again ona 
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Control Card # 3 


Column 


1-10 


11-20 


21-30 


Notation 


ACDAT 


EN2 


OMCJ 


zone 7 to centre j basis. If a 0 is specified, 
only the first set of cards will be read. If, 
however, actual consumer survey and/or distance 
data are also to be read, a 1 must be entered. 


Definition 


An indicator that, when assigned a value of 1, 
specifies that actual distance data are to be 
used, i.e. 2 sets of zone cards are to be read. 
If ACDAT is equal to 0, then both actual 
consumer survey and distance data will be read, 
i.e., all 3 sets of zone cards are included. 
ACDAT must equal 0 if the lambda search process 
is desired. If neither actual distance nor 
consumer survey data are to be read, i.e.; no 
search is to be conducted, the space for this 
indicator may be left blank. In this case 

CTRI should be specified as 0. If actual 
distance data are to be used but no lambda 
search is desired both CTRI and ACDAT must be 
assigned a value of 1. 


A numerical value representing the number of 
paired values relating size of retail centre 

to profit percentage. This number cannot 
exceed the values associated with SL and FACT 

in the DIMENSION statement. If no profit or 
percentage profit calculations are to be made 

a value of 1 must be entered for this indicator. 


An indicator that, when denoted by the value 1, 
causes the calculation of a component of the 
denominator of R, the coefficient of correlation, 
to be calculated. This value is to be used if 
the lambda search process is desired and RFIND 
equals 0, If, however, the opposite holds true 
and RFIND equals 1, OMCJ must be specified as 0. 
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OMPE 


DATA 


ASC 


An indicator that, when assigned a value of 0, 
causes the calculation and printout of the 
profit at any centre j and the percentage 
profit factor for that .centre to be omitted. 
If this step is to be included, however, this 
indicator must be assigned a value of 1. If 
OMPE is not equal to 0, then EN2 must be a 
value greater than 1. 


An indicator that, when assigned a value of 1, 
ensures the termination of the program. If 
several data decks are to be run simultaneously 
DATA must equal 1 in the last deck and 0 in 
each preceding deck. If only a single deck is 
employed DATA must equal 1. 


An indicator that, when assigned a value of Ie 
causes actual sales to be calculated within 
the program on the basis of given sales per 
Square foot figures. If actual sales by 
centre are known, they may be read directly 
into the program. In this case, ASC must be 
set equal to 0. This indicator does not have 
to be specified if SITEN is assigned a value 
other than 0, since in this case actual sales 
are arrived at somewhat differently. 


Parameter Specification Card 


Column 


1-10 


Variable 


ALASP 


Definition 


A predetermined value specified for the 
parameter lambda. If no search procedure is 
utilized, the first set of zone cards in the 
data deck will initially identify this 
predetermined value as ALA]. In subsequent 
runs, however, it is conceivable that the 
testing of alternative parameter values will 
be desired. To obviate the time consuming 
task of repunching the zone cards each time 
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a new parameter value is to be run the 
indicator ALASP may be used. The program, 
when run without a search procedure, will 
substitute for lambda the value specified by 
ALASPeiTeAlAlmissoritginal ly specified as 1.0 
On each zone card. If no change in the 
parameter value is desired ALASP should be 
set equal to ALAl. If the parameter search 
process is initiated, ALASP should coincide 
with the lower limit of lambda selected. 


Profit Factor Card(s) 


Column Variable Definition 


1-10 SL(K) A numerical value representing the size or 
scale of retail operation related to a given 
DLO itt race Lon. 


11-20 FACT(K) A numerical value representing the profit 
fraction associated with a given size or scale 
of retail operation. 


One such card is required for each paired combination of SL(K) 
and FACT(K). The total number of profit factor cards must thus 
correspond to the value previously given to EN2 on control card # 3. 
If EN2 is set at 1, then only 1 card is required with both SL(K) and 
PACT(K) specified as 0. 


Location Cards 


Column Variable Definition 
f-10 S(J) A numerical value corresponding to the size in 


square feet of a given retail centre j. 


11-20 XS(J) The X co-ordinate of a given retail centre j 
calculated from a grid map of the study area. 


21-30 YS(J) The Y co-ordinate of a given retail centre j 
calculated from a grid map of the study area. 
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These X and Y co-ordinates are used for the 
calculation of straight-line distances in the 
program. They are essential if no actual 
distance data is available. Where actual 
distance data is available and no straight-line 
output is desired, these co-ordinate values 
need not be specified. 


31-40 DS(J) A numerical value corresponding to the size 
increment. Chate 5 stombesapp lied toma given 
retail centre j if called for by the indicator 
SITEN On controle cardsye2e i tano incrementing 
is to be done this space must be left blank. 
That is, if SITEN equals 0, no value is to be 
entered for DS(J). 


41-50 SM(J) A numerical value representing the number of 
times the size increment DS(J) is to be applied 
to a given retail centre j if called for by 
the indicator SITEN. Again, SITEN cannot 
equal 0 if a value is to be specified. 


There will be one location card for each retail centre included 
in the study area. That is, the number of location cards must 
correspond to the value entered for ENJ on control card #1. It is 
helpful to arrange these cards in descending order of size of retail 
footage. This allows the program to be run more readily for various 
ranges of *centre sizes, 


Combined Sales Per Square Foot Cards 


Column Variable Definition 
1-5 — FS ES(0) An average retail sales per square foot value, 


i.e., FOOD and DSTM combined, associated with 
each centre in the study area. This value is 
to be specified only if site incrementing is 
desired. That is, SITEN cannot be 0. In this 
case, no data for AMKB(J), FOOD(J), DSTM(J), 
SLSF(J) or SLSD(J) can be included. A 
separate card is required for each centre. 
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mone Cards: Set # 1 


Column Variable Definition 


1-10 CT A numerical value corresponding to the total 
number of consumers or households in a given 
ZONemais. 


11-20 ALA] A numerical value representing the lower limit 
of the parameter lambda for a given zone i if 
a search procedure is being utilized. If no 
search is to be conducted, a predetermined 
value for lambda must be entered. It is 
simplest to assign ALAl a value of 1.0 as 
lambda can then be easily respecified by the 
indicator ALASP on the Parameter Specification 
Card. 


21-30 ALA2 A numerical value representing the upper limit 
of the parameter lambda for a given zone i if 
a search procedure is being utilized. If no 
search is to be conducted the space for this 
term must be left blank. 


31-40 BF A conversion factor on CT, the number of 
consumers or households in a zone i, representing 
retail expenditures for various products 
Originating in that zone. In effect this 
represents the average annual income per 
consumer or household available for retail 
spending in a zone i. If no sales distribution 
(by centre) is desired or no income data are 
available, a O should be noted. This will 
result in a printout of the customer 
distmipucvoneaonly. ein virtuallyvall practical 
applications of the model, however, it is most 
important to specify income data in some form. 


41-50 XN The X co-ordinate of the centre of population 
gravity of a given zone i, usually the 
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Zone Cards: Set 4-2 


Column 
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Variable 


C(J) 


estimated centroid of the zone, calculated 
from a grid map of the study area. 


The Y co-ordinate of the centre of population 
gravity of a given zone i calculated from a 
grid map of the study area. These X and Y 
co-ordinates are also critical for the 
calculation of straight-line distances in the 
program. They are essential if no actual 
distance data are available but may be omitted 
if straight-line output is not required. 


Definition 


A numerical value representing the actual 
number of shopping trips made from a given 
zone i to a given centre j. Provided that 
such data are available, each zone can have 
One corresponding consumer survey card(s). 
Each such card(s) represents the actual 
distribution of consumers from a given zone i 
to each centre j. A separate value for each 
centre j is required, with a maximum of 6 
columns allowed per entry. Data for as many 
as 12 centres can be accommodated on a single 
card, with remaining centres being entered on 
additional cards as necessary. The indicators 
ACDAT and CTRI must equal O and 1 respectively 
if such data are to be read. If no figures 
for actual consumer distribution are available, 
no card(s) is required. In this case, the 
space for ACDAT must be left blank if actual 
distance data are also not to be used, or 
conversely, assigned a value of 1 if such 

data are available. Correspondingly, CTRI 
should be either O or 1. 


Zone Cards: Set # 3 


Column Variable Definition 


1-72 TCH A numerical value representing the actual 
diatancesin trave P timeyscosty mites, etc. . 
from a given zone i to a given centre j. 


Provided that actual data are to be used, each 


zone must have one corresponding distance 
card(s). Each such card(s) specifies the 
actual distances from a given zone i to each 


centre j. A separate distance value is entered 


for each centre, with a maximum of 6 columns 
allowed per entry. Up to 12 distance values 


can be noted on a single card, with additional 


values entered on extra cards as necessary. 


Again, the indicator ACDAT on control card # 3 


must equal 0 if both actual consumer survey and 


actual distance data are to be used. If only 


actual distance figures are available, however, 


ACDAT must equal J]. If neither the second nor 


the third set of zone cards are to be used, 
the space for ACDAT should be left blank and 
CTRI should equal O. 


The total number of zone cards in each of the three possible 
Sets must correspond to the value specified for ENI on control 
card # 1. If actual consumer survey and distance data are to be 
utilized the following procedure must be adopted. Each zone card 
must be followed, first, by its corresponding consumer survey data 
(C(J)) card, and second by its corresponding distance data (T(J)) 
card. If only actual distances are to be employed, a distance 
data card must follow immediately after each zone card. The cards 
must be alternated for each zone in this manner. 


Actual Retail Sales Card(s) 


Column Variable Definition 
1-70 AMKB(J) A numerical value corresponding to the total 


annual dollar sales actually generated at a 
given retail centre j. Seven such values can 


OG N 


Floorspace Card(s): 


Floorspace Card(s): 


be accommodated per card with a maximum of 10 
columns allowed per entry. Additional values 
are to be specified on extra cards if necessary 
These values must be ordered so as to correspond 
to the array of location cards. The indicator 
ASC on control card # 3 must equal O if actual 
sales are given. Actual sales figures are 
usually not available, however, and must be 
estimated. If this is the case, this card(s) 

is to be omitted and the indicator ASC must 
equal 1. If SITEN is a value other than 0, 

this card(s) must also be omitted. 


Set # ] 


Definition 


A numerical value corresponding to the total 
retail floorspace in square feet at any given 
centre j devoted to the sale of food items. 


Set # 2 


Definition 


A numerical value corresponding to the total 
retail floorspace in square feet at any given 
centre j devoted to the sale of department 
store type merchandise. Values noted in each 
of the 2 sets of floorspace cards must 
correspond to the sequence of centres specified 
in the location cards. In each case seven 
values per card may be entered, with a maximum 
of 10 columns allowed per entry. Additional 
values must be specified on extra cards as 
necessary. No card(s) for either FOOD(J) or 
DSTM(J) can be included if SITEN has a value 
Other than 0. 
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Sales Per Square Foot Cards 


Column Variable Definition 
1-10 SESECd) A numerical value specifying the retail food 


sales per square foot assigned to each centre 
in the study area. 


11-20 SLSD(J) A numerical value specifying the retail DSTM 
Sales per square foot assigned to each centre 
in the study area. 


These variables provide one method of estimating actual annual 
Betail sales. If no distinction in sales per Square foot is to be 
drawn between any of the centres in the study area, identical 
values should be noted for each case, with a separate card included 
for each centre. However, it may be possible in some cases to 
determine individual sales per square foot values for each retail 
centre. 


One example of a possible basis of differentiation would be the 
criterion of planned versus unplanned centres. In this case, one 
card per centre is entered as before, but values need not be the 
Meer lf G€ither the indicator ASC is specified as 0, or SITEN is 
not equal to 0, no sales per square foot card(s) should be included 
with the input deck. If actual sales are to be estimated on the 
basis of these two variables, no input for AMKB(J) is necessary. It 
Should also be noted that if AVIND is specified as 1, it is necessary 
Beerepeat the last 3 data cards in identical sequence, i.e., 
micorspace card(s): set #1, floorspace card(s): set # 2, and 
Sales per square foot card(s). 


Description of the Output Variables 


It will be readily apparent from the input data that a wide 
variety of options are available to the user. Although the 
Selection of options will be determined in large measure by the 
nature of the information available for a given study area, the 
number of alternatives available still allows the user considerable 
Flexibility. Consequently, the program can test the implications of 
a wide number of situations. 
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It follows, as well, that there will be a wide range of output 
variables associated with the program. The exact form of the 
printed results is dependent entirely upon the choice of control 
card indicators and the nature of the input data. A comprehensive 
listing has been prepared of all possible output variables produced 
by the program. The variables are defined in the sequence in which 
they appear in the body of the source program. 


The nature of the application being made and the type of analysis 
to be attempted will determine the value to the user of each output 
variable. Obviously, not all output variables will appear in the 
results each time. The options specified by the control cards will 


determine those that are to appear. For the most part the output 
variables that do appear are given descriptive titles in place of 
program variable names. The equivalent forms are indicated in those 


instances where titles are employed in place of program variable 


names. 


Output Variable 
Sek) 


FACT(K) 


AV 


SW 


Definition 


The size or scale of retail operation related 
tova given profitefractione 


The profit fraction associated with a given 
size of retail operation. Both of these 
variables appear under the heading PROFIT 
FACTOR TABLE in the printout. 


A fractional value equivalent to 1 divided by 
the number of retail centres. 


A value corresponding to a component of the 
denominator of R, the coefficient of correlation. 


The total number of actual or observed 
consumers or households going to all centres 
from a given zone i, i.e., the sum of C(J). 


These variables are associated with the 
Fibonacci search procedure included in the 


Huff program. For specific use, see 
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A Programmed Solutton for Estimating Retail 
Sales Potenttals (Huff and Blue 1964). 


Neighbourhood Printout (Outputs) 


Output Variable Definition 
oT The number of consumers or households in a 


given zone i. 


XN The X co-ordinate of the centroid of population 
dVavietyvaot anvazone 1. 


YN The Y co-ordinate of the centroid of population 
gravity of any zone i. 


ALA] The initial value of the lower limit assigned 
to the parameter lambda if the search option 
is utilized. If no search procedure is employed, 
ALAl corresponds to a predetermined lambda value. 


ALA2 The initial value of the upper limit assigned 
to the parameter lambda if the search option 
is utilized. If no search procedure is 
employed, the value O appears. 


ALA3 The actual value of the parameter lambda for a 
given zone i. This value may be derived by 
means of the search procedure in which case it 
represents an optimum value. It may also 
represent an average lambda value calculated 
over all zones. Finally, if no search 
procedure is desired, it may be a value 
specified by the user, in which case it 
corresponds to ALAl. 


R The coefficient of correlation associated with 
a given value of lambda. This value of R may 
correspond to either of the 2 possible values 
associated with ALA3 if a search procedure is 
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BF 


EC(J) 


ECB(J) 


utilized. The indicator RFIND must be 
specified as 0 or R will not be calculated, 
appearing as 0 in the printout. R cannot be 
calculated if the search process is omitted as 
the calculation of this coefficient is 
dependent upon the input of actual consumer 
survey data. 


The average annual consumer or household income 
available for retail spending in a zone i. 


An integer number used to identify each retail 
centre j. This value appears under the heading 
SITE in the neighbourhood printout. 


The distance from a given zone i to a retail 
centre j. This value can represent either 
actual distance in terms of travel time, cost, 
etc., Or calculated straight-line distance. 
depending on the option selected and the data 
available. This variable is identified as 
DIST in the neighbourhood printout. 


The expected number of consumers or households 
originating in a zone i and shopping at a 
centre Jj. thts variable is identified as, CUcm 
in the neighbourhood printout. 


The expected retail spending originating in a 
zone i and terminating at a centre j. In 
effect, this value represents the expected 
number of consumers from zone i shopping at 
centre j multiplied by the retail spending 
available at zone i. This variable appears 
under the heading EXPEND in the neighbourhood 
printout. 


The total retail floorspace in square feet at 
a given centre j. This variable is identified 
as SIZE in the neighbourhood printout. 


eed ae 


XS(J) The X co-ordinate of a given centre j. This 
value appears under the heading XS in the 
neighbourhood printout. 


YS(J) The Y co-ordinate of a given centre j. This 
value appears under the heading YS in the 
neighbourhood printout. 


ET(Jd) The probability of a consumer or household in 
a given zone i shopping at a given retail 
centre Jj. "This variable is identified as 
HSHARE in the neighbourhood printout. 


ECT The sum of the EC(J) values for a given zone i. 
ECT thus represents the total number of 
consumers or households in a given zone and 
should correspond to the value for CT. 


ECBT The sum of the ECB(J) values for a given 
zone i. ECBT thus represents the total retail 
spending available in a given zone. 


ETT The sum of the ET(J) values for a given zone i. 
The sum of these individual probabilities must 
dlWaysebewiln0n his: 15 the last variable to 
appear in the neighbourhood printout. 


AVLAM An average value of lambda calculated by 
summing the optimum parameter values for each 
zone and dividing by the number of zones. 


PROF (J) The dollar profit associated with a given 
retail centre j. This variable is identified 
CSerRO Bl  Aaneechicap ri mtout. 


PERC (J) The percentage profit factor, determined from 
the profit factor table, associated with a 
given centre j. This variable is identified 
as PERCENT ine the printout. 
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Summary Site Printout (Outputs) 


Output Variable 
J 


Si) 


FOOD(J) 


DSTM(J) 


XS(J) 


YS(J) 


TMKT (J) 


Definition 


An integer number used to identify each retail 
centre j. This value appears under the heading 
SITE in the summary site printout. 


The total retail floorspace in square feet at 
a given centre j. This variable appears under 
the heading SIZE in the summary site printout. 
The sum of these values is denoted in the 
program as Tl. 


The floorspace in square feet devoted to the 
sale of food itemseateay given centre yj. sinha 
variable is identified as FOOD in the summary 
site printout. The sum of these values is 
denoted by T2. 


The floorspace in square feet devoted to the 
sale of department store type merchandise at 

a given centre j. This variable is identified 
as DSTM in the summary site printout. The sum 
of these values is denoted by T3. 


The X co-ordinate of a given centre j. This 
value appears under the heading XS in the 
summary site printout. 


The Y co-ordinate of a given centre j. This 
value appears under the heading YS in the 
Summary site printout. 


The total number of consumers or households 
from all zones shopping at a given retail 
centre j. This variable appears as CUST in 
the summary site printout. The total number 
of consumers or households shopping at all 
centres is designated as T4. 
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TMKB (J) 


AMKB(J) 


RATIO(J) 


ACTRAT(J) 


FRTMK(J) 


SLSMK(J) 


The total retail sales generated by a given 
centre j as calculated by the model. This 
Varivable is identified as CAL SLS dn the 
Summary site printout. The sum of the 
calculated retail sales generated by all 
centres in the study area is designated as T5. 


The actual sales generated by a given centre j. 
This variable is identified as ACT SLS in the 
summary site printout. The sum of the actual 
retail sales generated by all centres is denoted 
DVO. 


The ratio of calculated sales to total 
floorspace at a given centre j. In effect, 
this variable represents the sales per square 
foot generated at a centre j as calculated by 
the model. This variable is identified as 
CS/FT in the summary site printout. The 
average calculated sales per square foot for 
all centres is denoted in the program by T7. 


The ratio of actual sales to total floorspace 
at a given centre j. This variable represents 
the sales per square foot actually generated 
at a centre j and is identified as AS/FT in 
the summary site printout. The average actual 
sales per square foot for all centres is 
identified as T8. 


The percentage of total consumers or households 
in the study area patronizing a given centre j. 
This variable is identified as HH SH in the 
Summary site printout. The sum of these 
percentages is denoted by T9. 


The percentage of total calculated sales 
generated by a given centre j. This variable 
is identified as CA$SH in the summary site 
printout. The sum of these percentages is 
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designated as T10. 


ASLSMK(J) The percentage of total actual sales generated 
by a given centre j- This variable, is 
identified as AC$SH in the summary site output. 
The sum of these percentages is identified in 
the programas Vll. This is the last variable 
to appear in the summary site output table. 


THE CULLEN PROGRAM 


The second program is based on the Cullen formulation of the 
retail gravity model incorporating the Wilsonian concepts of entropy. 
The version presented here is actually an extensive revision of a 
program developed by Ian Cullen at University College, London. The 
sequence of presentation is identical to that used for the presentation 
of the Huff version in the first part of this appendix. 


The Source Program 


A copy of the source program will be found in Annexure B. 
Explanatory notes and a flow diagram have been added in the same way 
as in the Huff version in Annexure A. This program is also written 
in FORTRAN IV. 


The notational form of the Cullen version is 


In the program version it appears as 


ALPHA (A(J)) - LAMBDA (D(I,J)) 


Pie | e 
S(I,J) = ALPHA (A(J)) - LAMBDA (D(I,J 


The translation of the notational form of variables to the 
program form is explained in Fig. A-2. 
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NOTATIONAL AND PROGRAM VARIABLES - CULLEN 


Definition 


The total retail expenditure available 


Mathematical 
Notation 


Program 
Notation 


lied ON Cale 


The sales attracted to a centre JOfrom 


azZOne 1. 


A measure of the attractiveness ot a 


Centreoj; in this case total annual 


retail sales. 


A measure of travel distance, time or 


cost from a zone i to a centre ae 


The number of retail centres in a 


given study area. 


A parameter estimated empirically to 


reflect the attractiveness of a 


Cencres je 


A parameter estimated empirically to 
werlect the effect 01; travel distance, 


time or cost on various kinds of 


Shopping trips. 


Fig. A-2: Notational and Program Variables - Cullen 


Sf 


a 


a 


Description of the Input Data 


The instructions for the preparation of the input data in 
punch card form are listed below. The cards must be arranged in 
the same sequence in which they are described. The modified Cullen 
program provides the user with a variety of procedural options, each 
of which is fully outlined as it appears. 


Study Area Card 


Column Variable Definition 
1-3 N An integer value representing the number of 


residential zones in the study area. 


4-6 M An integer value representing the number of 
retail centres in the study area. 


Control Card # 1 


Column Notation Definition 
1-2 NA An indicator that, when assigned a value of 1, 


specifies that a calibration run is to be made. 
If a calibration run is not desired, column 2 
must be left blank. 


3-4 NB An indicator that, when assigned a value of 1, 
specifies that a study year simulation run is 
to be made. If such a run is not desired, 
column 4 must be left blank. 


5-6 NC An indicator that, when assigned a value of 1, 
specifies that a projection run is to be made. 
If such a run is not desired, column 6 must 
be left blank. 


7-8 NE An indicator that specifies the number of 
periods into which a given time span is to be 
subdivided if a projection run is to be made. 
For example, if the base year was 1971 and 
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17-18 


19-20 


NF 


NG 


the design year 1976, a value of 5 would 
result in a printout for each year in this 
time period, i.e., 1972-1976. If NA or NB 
areuespecifivedsas) lowthismespace snouldebemle ht 
blank. 


A numerical value identified as the constant 
of retail consumption. This is a fixed value 
used to convert consumer or household income 
to its retail spending equivalent. If left 
blank or specified as 0, this constant is 
calculated automatically from the input data. 
However, more reliable results are likely to 
be obtained if an actual value is provided. 
If consumer or household retail spending data 
are available the constant must be specified 
as 1.0. If either of the options associated 
with NG or NH are to be adopted C must also 
be specified as 1.0. 


An indicator that must be specified as 1 if 
travel» cost=datesane sinnuteas straight-line 
centimetre distances. If this is not the 
case, this column must be left blank and some 
true measure of travel cost used for D(I,J). 


An indicator that, when assigned a value of 1, 
specifies that base year sales figures for M-1 
of M centres are to be adjusted, each in 
proportion to the overall increment in total 
consumption implied by design year population 
and income data. This option should be 
specified only for a study year simulation run 
in which the last, i.e., Mth, sales figure 
input is a design year estimate for a proposed 
centre and all other sales figures relate to 
the base year. If this option is not desired 
NG should be left blank. 


=O 


ei 22 


NH 


Control Card # 2 


Column 


1-3 


4-6 


Notation 


ASC 


DATA 


Trip Length Card 


Column 


1-8 


Notation 


AST 


An indicator that, when assigned a value of 0 
or left blank, specifies that base year sales 
figures are to be adjusted as above FO.Veugals ling 
lve. MNecentres. pif NHois specified as 0, NG 
should equal 1. This option is to be used if 
a study year simulation run is being made on 
the basis of base year sales figures to test 
Eiemet jects otenoteintroducingsa Proposed centre. 
If this option is not desired NH should be 
left blank. Summarizing, if neither option is 
desired both NG and NH must be left blank. If 
NG is assigned a value of 1 and NH is blank or 
0, all centres will be adjusted. If both 
indicators are specified as 1, M-1 centres 
will be adjusted. 


Definition 


An indicator that, when assigned a value of ils 
causes actual sales to be calculated within 
the program on the basis of given sales per 
Square foot figures. If actual sales by 
centre, i.e., A(J), are known they may be 
readudirectly into the program. In this case, 
ASC must be specified as 0. 


An indicator that, when specified as 1, ensures 
the termination of the program. If a number 

of data decks are to be run simultaneously, 
DATA must equal 1 in the last deck and 0 in 
each preceding deck. If one deck only is used, 
DATA must be specified as 1. 


Definition 


A numerical value representative of the 
observed average shopping trip length in the 
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study area. If no such value is desired, a 
blankivcard mus Gabes sertedsat this pone. 


Travel Speed Card 


Column Notation Definition 


1-8 V A numerical value representing the average 
travel speed in miles per hour over the study 
area. A figure must be specified unless some 
actual measure of travel cost is utilized for 
D(I,J), in which case a blank card must be 
inserted. If Vis left blank, the sindicator 
NF on control card # 1 must also be left blank. 


Distance Cards 


Column Variable Definition 


1-72 DCT sd) A numerical value representing the cost of 
travelling from the centroid of a given 
residential zone i to a retail centre j. This 
value may simply be a straight-line centimetre 
distance measured from a standard grid map of 
the study area. In this case a distance 
conversion procedure may be applied by setting 
the indicator NF equal to 1. Alternatively, 
actual travel costs may be specified in terms 
of dollars, miles, time. etce. snouldesucn 
data be available. In this case, both NF and 
V should be left blank. As in the modified 
Huff and Blue program, a separate value is 
noted for each centre, with a maximum of 6 
columns allowed per entry. This allows a 
maximum of 12 such values to be specified per 
card, with additional values noted on 
successive cards if needed. The cards must 
be arranged consecutively by zone with the 
total number of cards corresponding to the 
value of N specified in the study area card. 
The keypunching of distance data is one of 
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Population Card(s) 


Column Variable 
1-80 P(I) 


Floorspace Card(s): 


Column Variable 
1-70 FS(J) 


Floorspace Card(s): 


Column Variable 
1-70 FOOD(J) 


the more tedious and time-consuming tasks 
associated with the preparation of the input 
deck. It should be noted that the format used 
allows the same distance cards to be used in 
either of the 2 programs. Care must be taken, 
however, to (i) properly Specify the associated 
control indicators, and (ii) to arrange the 


cards correctly. 


Definition 


A numerical value representing the number of 
consumers or households in a given residential 
zone i. Separate values must be noted for 

each zone, with a maximum of 10 columns allowed 
per entry. A maximum of 8 zonal population 
figures can thus be specified per card, with 
values being punched consecutively by zone 
If either a 


calibration run or a study year simulation 


On aS many cards as required. 


run is being made these values correspond to 
study year zonal population. If, however, a 
projection run is being made they must 
correspond to design year population figures. 
In this latter case, base year data are 
included at the end of the data deck. 


Set # ] 


Definition 


A numerical value corresponding to the total 
retail floorspace in square feet at any given 
Cen Urea. 


Seu. 7 ac 


Definition 


A numerical value corresponding to the total 


Seetis 2 


retail floorspace in square feet at any given 
centre j devoted to the sale of food items. 


Floorspace Card(s): Set # 3 


Column Variable Definition 


1-70 DSTM(J) A numerical value corresponding to the total 
retail floorspace in square feet at any given 
centre j devoted to the sale of department 
store type merchandise. 


Each set of floorspace cards is punched under an identical 
format. In each case up to seven values may be entered per card, 
with 10 columns available for each entry. Additional values are to 
be entered on extra cards as required. Values must be entered 
consecutively by centre. Since the format for the variables FOOD(J) 
and DSTM(J) coincides with that employed in the modified Huff and 
Blue program these data cards can be substituted in either of the 
two programs. 


Sales Per Square Foot Cards 


Column Variable Definition 


1-10 SESE (uy) A numerical value specifying the retail food 
sales per square foot assigned to each centre 
in the study area. 


11-20 SLSD(J) A numerical value specifying the DSTM sales 
| per square foot assigned to each centre in 
the study area. If actual sales are to be 

estimated on the basis of these 2 variables, 

no input for the variable A(J) is necessary. 


This input allows actual retail sales to be estimated for each 
centre on the basis of sales per square foot data. If this option 
is to be employed. one card must be punched for each centre. If no 
distinction is to be made in sales per square foot, identical values 
should be entered on each card. Differential values may be used if 
individual sales per square foot estimates are available for each 
centre. These sales per square foot cards may also be substituted 


oe A 


Income Card(s) 


Column 


Actual 


Column 


1-70 


Variable 


into either of the two programs, as identical 
formats are required. If, however, the 
indicator ASC is specified as 0, no sales per 
square foot cards are required. 


Definition 


A numerical value representing the average 
annual household or consumer income in dollars 
for a given residential zone i. This value is 
reduced to its retail spending equivalent by 
the consumption constant C entered on control 
card #1. If, however, actual income budgeted 
for retail expenditure by zone is known, it 
may be specified accordingly. In this case, 

C must be specified as 1. If either a 
calibration run or a study year simulation 

run is being made, these data correspond to 
study year income figures. If, however, a 
projection run is to be made, these data must 
correspond to design year income. In this 
case, base year income data are included at 
the end of the data deck. 


Retail Sales Card(s) 


Variable 


A(J) 


Definition 


A numerical value corresponding to the total 
annual dollar sales actually generated at a 
retail centre j. Seven such values can be 
entered per card with a maximum of 10 columns 
available for each entry. Data must be 
entered consecutively by centre and additional 
values must be specified on extra cards as 
necessary. If such data are available, ASC 
should be assigned a value of 0. If, however, 
actual sales for each centre are to be 
calculated within the program, ASC should 


=a 8)i1' 5: = 


equal l.5 Tnethis vcasesmno sactualesales 
card(s) is required. This data card(s) may 
also be substituted into either of the two 
programs used in this study. 


Calibration Parameter Card 


This card is required only if the indicator NA on control 
card # 1 has a value of 1. If no value for NA is specified then 


this card is to be omitted. 


Column 


1-9 


10-18 


19-20 


21-29 


30-38 


39-40 


Variable 


AL] 


ALI 


TAL 


ALAMI 


ALAMI 


ILAM 


Definition 


The initial value assigned to the attraction 
parameter, alpha, in the calibration run. 


A numerical value specifying the size of the 
alpha increment selected for the calibration 
run. 


An integer value specifying the number of 
increments to be used for the attraction 
parameter, alpha, in the calibration run. 


The initial value assigned to the distance 
parameter, lambda, in the calibration run. 


A numerical value specifying the size of the 
lambda increment selected for the calibration 
run. 


An integer value specifying the number of 
increments to be used for the distance 
parameter, lambda, in the calibration run. 


Parameter Specification Card 


This card is required only if either of the indicators NB or 


NC on control card # 1 is assigned a value of 1. This card is not 


required if a calibration run is to be made. 
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Column Variable Definition 


1-9 ALPHA The value assigned to the attraction parameter, 
alpha, if either a study year simulation run 
or a projection run is being made. 


10-18 LAMBDA The value assigned to the distance parameter, 
lambda, if either a study year simulation run 
or a projection run is being made. These 
parameter values are normally derived by means 
of a calibration run using base year data. 


Base Year Population Card(s) 


Column Variable Definition 


1-80 PB(I) A numerical value corresponding to the number 
of consumers or households in a given 
residential zone i in the base year. The 
format used is identical to that for the 
Variable P(I):. This card(s) is required only 
if a projection run is to be made. 


Base Year Income Card(s) 


Column Variable Definition 


1-80 YB(1) A numerical value corresponding to the average 
annual consumer or household income in dollars 
for a residential zone i in the base year. 
Again, this value may be reduced to its retail 
spending component by the consumption constant 
Cee lfeactuals retail) spending data are 
available, C must be assigned a value of 1. 
The format used is identical to that for the 
variable Y(I). This card(s) is required only 


itmoaprouection. run 81s) to be made. 


Description of the Output Variables 


The wide range of options available in the modified Cullen 
program is associated with an equally wide range of output variables. 
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The comments made about the output variables in the modified Huff 
and Blue version apply equally to the Cullen outputs (Pages 101-102). 
A complete listing of these variables follows. 


Calibration Run Outputs 


Output Variable 
C 


AST 


DSQ 


Definition 

The constant used to convert total household 

or consumer income to its retail spending 
equivalent. This may be either a predetermined 
value specified on control card #1, or it may 
be a value calculated within the program. This 
value is denoted as CONSUMPTION CONSTANT in the 
calibration run printout. 


The observed average shopping trip length in 
the study area. This is merely the value noted 
on the trip length card. This value is denoted 
as INPUT TRIP LENGTH in the calibration run 
printout. 


An integer value used to number the output 
printed with each iteration specified in the 
calibration run. 


The sum of the squares of the differences 
between actual and calculated sales for each 
centre j. This value is calculated in 
Subroutine DIFFSQ. It appears under the 
heading DIFFSQ in the calibration run printout. 
Theoretically, the parameter values associated 
with the minimum DSQ value best approximate 
reality. 


The value of the coefficient of correlation 
calculated with each iteration specified in 

the calibration run. This value is calculated 
in subroutine CORR and is intended to reflect 
the degree of correlation between calculated 
and actual sales. It appears under the heading 


= e.o ee 


CORR in the calibration run printout. 


ASTL An average shopping trip length calculated 
for the entire study area with each iteration 
specified in the calibration run. This 
variable appears under the heading TRIP LENGTH 
in the calibration run printout. 


ALPHA The value of the attraction parameter, alpha. 
A separate value is printed for each iteration 
specified in the calibration procedure. These 
values appear under the heading ALPHA in the 
calibration run printout. 


LAMBDA The value of the distance parameter, lambda. 
Again, a separate value is printed for each 
iteration specified in the calibration 
procedure. These values appear under the 
heading LAMBDA in the calibration run printout. 


J An integer number used to identify each retail 
centre j. This output is denoted by the heading 
SITE in the calibration run printout. 


A(J) The actual sales generated by a given retail 
centre j. This variable may be either 
estimated within the program or read directly 
from the data deck. It is identified by the 
heading INPUT SALES in the calibration run 
printout, 


T(J) The total retail sales generated by a given 
centre as calculated by the program using the 
parameter values associated with each iteration 
specified in the calibration run. This variable 
appears under the heading SYNTHESIZED SALES in 
the calibration run printout. 


= Sato = 


Study Year Simulation Run Outputs 


Output Variable 
DSQ 


FS(J) 


FOOD(J) 


DSTM(J) 


SUMT (J) 


Definition 


The sum of the squares of the differences 
between actual and calculated sales for each 
centre). This Valuew1s, calculatedsin 
subroutine DIFFSQ and provides an indication 
of the goodness of fit obtained by the model. 


An integer number used to identify each retail 
centre j. This output is denoted under the 
heading SITE in the summary site printout. 


The total retail floorspace in square feet at 
a given centre j. This variable is denoted as 
SIZE in the summary site printout. The sum of 
the total floorspace figures is designated in 
the program as Tl. 


The floorspace in square feet devoted to the 
sale of food items at any centre j. This 
variable appears under the heading FOOD in the 
Summary site printout. The sum of the food 
floorspace figures is designated as T2. 


The floorspace in square feet devoted to the 
sale of department store type merchandise at 
any centre j. This variable appears under the 
heading DSTM in the summary site printout. 

The sum of the DSTM floorspace figures is 
designated in the program as T3. 


A value representing the average calculated 
trip length from all zones to a given centre j. 
This value reflects trip distances weighted by 
their associated customer flows. The average 
trip length over the entire study area is 
designated as T12 while the average trip 
length to each centre appears under the 
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TMKT (J) 


RATIO(J) 


ACTRAT(J) 


heading AVTL in the summary site printout. 


The total number of consumers or households 
from all zones shopping at a given centre j. 
This value is denoted by CUST in the summary 
site printout. The sum of these values is 
designated as T4 and should correspond to the 
study area population in terms of either 
consumers or households. 


The retail sales generated by a given centre j 
as calculated by the program. This value is 
referred to as synthesized sales and appears 
under the heading SYN SALES in the summary 

site printout. The sum of these values, that 
is, the total sales generated by all centres 

in the study area as calculated by the program, 
isenocvedsas lo. 


The actual sales generated by a given retail 
centre j. This variable is denoted as ACT 
SALES in the summary site printout and may be 
either estimated within the program or read 
directly from the data deck. The total actual 
sales generated by all centres is noted in 

the program as T6. 


The ratio of synthesized or calculated sales 
fomtotalp¢loorspace ate agiven centre j. ~ihis 
variable thus represents the sales per square 
foot generated at a centre j as calculated by 
the model. This variable is denoted as CS/FT 
in the summary site printout. The average 
calculated sales per square foot for all 
centres in the study area is designated as T/7. 


The ratio of actual sales to total floorspace 


ateasdiven centre J. I[his variable represents 
the sales per square foot actually generated 


=) L2li = 


FRTMK(J) 


SLSMK(J) 


ASLSMK(J) 


YC(1)) 


ALPHA 


LAMBDA 


at any given centre and is denoted by AS/FT 
in the summary site printout. The average 
actual sales per square foot for all centres 
within the study area is denoted Dye ioe 


The percentage of total consumers or households 
in the study area shopping at a centre ties Pots 
variable is denoted by HH SH in the Summary 
site printout. The sum of these percentages 

is designated in the program as 19, 


The percentage of total calculated sales 
generated by any centre j. This variable is 
denoted by CA$SH in the summary site printout. 
The sum of these percentages is noted in the 
program as T10. 


The percentage of total actual sales generated 
by a given centre j. This variable appears as 
AC$SH in the Summary site printout. The sum 
of these percentages is denoted by Tll. 


An integer number used to identify each 
residential zone i in the listing of zonal 
retail spending. 


The total amount budgeted for retail spending 

by all consumers or households in a residential 
zone i. This variable appears under the heading 
ZONAL RETAIL SPENDING. 


The value of the consumption constant used to 
convert total household or consumer income to 


its retail spending equivalent. 


The value of the attraction parameter, alpha, 
entered on the parameter specification card. 


The value of the distance parameter, lambda, 


ma 4a 


entered on the parameter specification card. 


osu ) The retail spending from a given residential 
zone i allocated to a centre j. These values 
are printed in a matrix format under the 
heading CASH FLOWS TO. 


Sol, ) The number of consumers or households 
travelling from a given residential zone i to 
shop at a given centre j. These values also 
appear in a matrix format and are identified 
by the heading CUSTOMER FLOWS TO. 


Projection Runs 


The printed output from this version of the program is not 
entirely suitable for making projection runs. Projection runs are 
intended to allow the impact of new centres to be assessed at 
different points in time. That is, the existence of a new centre 
may be specified for any given year within a defined projection 
period. Consequently, it becomes possible to test numerous 
alternatives within a temporal framework. Certain additions to the 
body of the program are necessary, however, if such a run is to be 
made. One simple illustration is discussed below. The new cards 
that will be required are denoted by asterisks and appropriate 
explanations are provided throughout. 


An Illustration 


Assume that a projection is to be made with base year population 
and income data for 1971. Assume also that NE is specified as 5, 
that there are 6 shopping centres in the base year, and that a new 
centre is to be added in the third year of the projection period. 
In this case output would be provided for each of five years, 1972 
through to 1976, with a new centre being introduced in 1974. The 
following additions would have to be made to the source program. 


140 CONTINUE 123 - 


Ker E IL 20 
Le (ke GEQNE GO m1.0 7190 


Tit Keir 0) N26 


Alba Ks Gees. ON, 


The value assigned to K specifies the 

year in which the new centre will be 
added. In this Cases ih Ketse less than 
3.0, the original number of centres will 
be retained. There will thus be 6 centres 
specified for 9/72 and 1973. 


Once K is equal to 3.0 or greater, the 
new centre is introduced. The number of 
centres, M, has been increased by 1 for 
the period 1974-1976. 


*IF (K.GE:3310)FS(7)=200000 


Since certain information about the new 
centre cannot be provided in the data 
deck, it must be presented at this 
point. Fi Gsitwevo besspecifieduaismtne 
total floorspace in square feet of the 
new centre. 


FIP (KeGE3 40) ROOD(7 = 30000 


The floorspace devoted to the sale of 
food items must also be specified for 
the new centre. 


PLP CK GE 30) DST Mira) len OU0 


Similarly, the floorspace devoted to the 
sale of department store type merchandise 
is also indicated for the new centre. 


*IF (K.GE.3.0)A(7)=13500000 


Finaldy.theractual sales likely to be 
generated by the new centre is specified. 
This will normally be calculated by 
applying estimated sales per square foot 


eeai240)- 


to the floorspace figures anticipated. 


DO 142 I=1,N 
MA=1 
*MB=6 The dimensions of the cash flow matrix 
are specified for the original number 
of centres. 
*IF (K.GE.3.0)MB=7 inewmatrix ts adjusted to include the 


new centre at the time of its introduction. 


172 CONTINUE 


MA=1 
*MB=6 The dimensions of the customer flow 
matrix are similarly specified. 
mrm kK. GE.3.0)MB=7 The customer flow matrix is adjusted to 


include the new centre. 


176 CONTINUE 


The value entered for M on the study area card must correspond 
to the number of retail centres in the design or target year. 
Distance data must also be punched for all centres being considered 
in the project run. Similarly, one sales per square foot card 


ah 5) = 


must be included for every centre. Finally, floorspace figures for 
any centres other than those found in the base year must be entered 
as 0. his applies to each set of floorspace cards en: ESi(d ie 
FOOD(J) and DSTM(J). The output variables associated with the 
projection run are identical to those described for the study year 
Simulation run, 


The projection run affords the user with considerable 
flexibility in testing alternatives. Although this Simple example 
has illustrated the addition of only one new centre, it is possible 
to introduce any number of proposed centres in any desired time 
sequence if the appropriate cards are added to the body of the 
program. These cards must be removed, however, if a calibration or 
a study year simulation run is to be made. 
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